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ABSTRACT 

The  following  papers  are  presented  in  this  report: 
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I.  INTRODUCTION 


A.  OBJECTIVE 

The  aim  of  the  work  reported  here  is  to: 

1.  Review  the  available  theoretical  and  experimental  data  on  the 
reactions  that  control  the  electron  density  in  missile  wakes 
and  in  the  environment  surrounding  nuclear  bursts. 

2.  Supplement  this  by  calculating  the  cross  sections  (reaction 
rates)  for  specific  reactions  of  importance  in  these  problem 
areas.  (In  its  initial  stages,  the  program  has  studied  re¬ 
actions  in  atmospheric  gases  such  as  dissociative  recombina¬ 
tion  [Example:  N0+  +  e  -*■  N  +  0]  and  related  processes.) 

3.  Compute  other  reactions  as  requested  by  ARPA. 

4.  Develop  general  computational  techniques  designed  to  improve 
the  accuracy  and  reliability  of  calculations  of  this  type. 

B.  RELEVANCE  TO  THE  DEFENDER  PROGRAM 

The  program  at  Defense  Research  Corporation  is  directed  towards 
the  analyses  of  reaction  rates  that  are  directly  applicable  to  several 
critical  problem  areas  in  the  Defender  program.  Some  of  these  are: 

1.  Radar  return  from  missile  wakes. 

2.  Nuclear  fireball  clutter  and  scintillation  effects. 

3.  The  expansion  of  nuclear  debris  at  high  altitudes. 

4.  Atmospheric  disturbances  introduced  by  nuclear  weapons. 

The  objective  at  Defense  Research  Corporation  has  been  to  devise 
practical  computational  techniques  for  estimating  the  cross  sections  of 
the  critical  reactions  connected  with  these  problems.  The  personnel 
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assigned  to  tills  task,  well  qualified  by  training  and  experience  to  under¬ 
take  this  work,  include: 

Dr.  B.  A.  Lippmann,  Principal  Investigator  (half-time). 

Dr.  T.  F.  O'Malley  (full  time). 

Dr.  P.  J.  Redmond  (half-time). 

C.  SYNOPSIS  OF  WORK  ACCOMPLISHED 

At  the  present  time,  the  work  is  proceeding  along  two  lines.  The 
main  emphasis  is  being  placed  on  numerical  analysis  of  specific  reactions. 
However,  auxiliary  formal  studies  relating  to  questions  arising  in  the 
course  of  attempting  the  numerical  calculations  are  also  being  carried  out. 

While  some  aspects  of  the  work  have  not  yet  been  completed,  or  are 
still  in  the  process  of  being  written  up,  a  fair  sampling  of  the  work  ac¬ 
complished  during  the  1966  period  is  contained  in  the  list  of  reports  and 
journal  articles  given  in  the  Table  of  Contents.  Copies  of  all  these 
documents  are  attached.  In  addition,  technical  assistance  in  connection 
with  the  ARPA  A&M  Physics  Program  has  been  furnisho-  as  requested. 
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II.  COMPANY  REPORTS  AND  JOURNAL  ARTICLES 


UNCLASSIFIED 


5 


UNCLASSIFIED 


THIS  PAGE  INTENTIONALLY  BLANK 


UNCLASSIFIED 


UNCLASSIFIED 


OPERATOR  FOR  TIME  DELAY  INDUCED 
BY  SCATTERING* 

by 

B.  A.  Lippraann 


ABSTRACI 

The  operator  that  gives  the  time  delay  induced  by  a 
scattering  process  is  exhibited  explicitly. 


It  is  well  knf'wn  that  the  time  delay  induced  by  a  scattering  process 

1  2 

is  given  by  the  energy  derivative  of  the  phase  shift.  Some  time  ago, 
we  observed  that  this  theorem  can  be  put  in  the  "suggestive"  form  S  *rS, 
where  r  is  interpreted  as  an  "operator,"  equivalent,  in  an  energy  repre¬ 
sentation,  to  differentiation  with  respect  to  the  energy.  It  is  the  pur¬ 
pose  of  this  paper  to  make  this  relationship  precise  by  exhibiting  the 
time-delay  operator  explicitly. 

ThtL  such  an  operator  might  exist  is  suggested  by  the  action  and 
angle  variables  of  classical  theory.  For,  if  the  action  variable  is 
identified  with  the  energy,  the  conjugate  angle  variable  has  the 


*  Published  in  Phy3.  Rev.  151,  1023  (1966). 

i 

This  appears  to  have  first  been  shown  quantura-mechanically  by  E.  P. 
Wigner,  Phys.  Rev.  98,  145  (1955),  although  relationships  of  this  sort 
are  well  known  in  filter  theory,  where  the  frequency  derivative  of  ':he 
phase  characteristic  gives  the  time  required  for  a  pulse  of  energy  to 
pass  through  the  filter  f  Radiation  Laboratory  Series  (McGraw-Hill  Book 
Company,  Inc.,  New  York,  1948),  Vol.  8,  p.  155].  A  similar  relation 
describes  the  spatial  displacement  of  a  wave  packet  upon  reflect ' jn  by 
a  plane  interface  whose  reflection  coefficient  is  a  function  of  the 
wave  number  parallel  to  the  interface,  L.  M.  Brekhovskikh,  Waves  in 
Layered  Media  (Academic  Press  Inc.,  New  York,  1960),  p.  105. 

^  See  F.  T.  Smith,  Phys.  Rev.  118,  349  (1960),  Eq.  (45). 
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equation  ot  motion 


3H 

3J 


1 


which  integrates  to 


♦  j(t) 


t  +  (0) 


(1) 

(2) 


One  would  expect  that  the  first  equation  would  have  as  its  quantum- 
mechanical  analog  the  conve.''*' ional  commutation  relation  that  holds  for 
conjugate  variables,  while  the  second  equation,  transcribed  to  quantum 
mechanics,  would  relate  one  of  the  quantum-mechanical  operators  to  the 
time  parameter.  It  is  shown  below  that  these  expectatior  .  can  be  realized. 


We  first  introduce  the  operator  for  the  magnitude  of  the  particle 
momentum 


P 


(3) 


and  the  operator  represer  .ing  the  projection  of  the  radial  coordinate 
upon  the  momentum 

rp  *  (1/P/(pt)  .  (4) 

3 

These  are  analogous  to  the  operators  introduced  by  Dirac  except 
that  the  roles  of  position  and  momentum  have  ’  sen  exchanged. 


As  in  Ref.  3,  we  easily  find  that  these  operators  satisfy  the 
commutation  relation 

[r  p]  -  i  (5) 

P  » 

and  that  the  operator  r^  +  i/p,  which  is  Hertnitian,  is  canonically  con¬ 
jugate  to  p  . 


P.  A.  M.  Dirac,  The  Principles  of  Quan^m  Mechanics  (Oxford  University 
Pr^ss,  New  York,  1958),  p.  152. 
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The  time-delay  operator  t  is  defined  by 


T  ■  !;(rp  +  i) +  ( rP  +  i • 


For  the  free  particle  Hamiltonian,  **  p  '2m,  we  find 


[t,Ho] 


i  . 


(6) 


(7) 


Although  t,  as  constructed  in  Eq.  (6)  above  is  formally  Hermitian, 
some  care  must  ba  *.aken  in  developing  its  properties,  because  of  the  fact 

4 

that  Hq  has  a  spectrum  that  is  limited  to  positive  values  (0,») ,  .  w- 

ever,  this  circumstance  does  not  impair  the  interpretation,  which  follows 
directly  from  Eq.  (7),  of  t  as  an  energy  derivative  in  a  representation 
in  which  Hq  is  diagonal:  r-»i3/3E  . 

Calling  the  operator  defined  in  Eq.  (6)  t(0),  we  can  also  introduce 
the  dime-dependent  operator 

iH  t  -iH  t 

T(t)  -  e  °  Tro)e  °  (8) 

-  t  +  t(0)  . 


The  relations  given  in  Eqs.  (7)  and  (8)  above  are  the  quantum  me¬ 
chanical  analogs  of  the  classical  Eqs.  (1)  and  (2). 

The  application  to  the  calculation  of  time  delays  is  made  by  con¬ 
sidering  the  matrix  element 

(^a(t),T(0)il>a(t)  i  ,  (9) 

where  the  state  vector  (t)  is  the  scattered  wave  packet,  evaluated  for 

fl 

times  after  the  scattering  process  has  been  completed.  The  time  depen¬ 
dence  of  this  state  vector  is  then  given  by  the  f ree-particle  Hamiltonian 
“iH  t 

alone:  i|i  (t)  ■  e  0  ^  (0):  here,  the  last  factor  represents  the 

&  3 


For  example,  the  operator  i  has  no  eigenfunctions;  see  W.  Pauli,  Hand- 
buch  der  Physlk  (Sprlnger-Verlag,  Berlin,  1958),  Vol.  5/1,  p.  63. 
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scattered-state  vector  extrapolated  back  to  zero  time. 

In  the  following,  both  the  incident  and  scattered  wave  packets  are 
assumed  to  be  normalized  to  unity.  Since  the  scattered  state  is  con¬ 
nected  with  the  initial  state  vector  *ft(0)  by  the  S  matrix,  (9)  can  be 
evaluated  as  below: 

,  iH  t  -iH  t 

(V°>»e  0  T(0)e  0  vo)) 

■  ('V0),|t  +  T<0)!  V0)) 

"  t  +  ( V'0),S-1T(0)S*a(0)).  (10) 

If  the  S  matrix  has  the  form  S  -  e216(E),  in  an  energy  representa¬ 
tion  the  time  delay  operator  becomes  equivalent  to  energy  differentiation, 
and  wc  find,  as  a  final  result 

(*a<t>,T<0>*a<t))  -  t  -  2  ^||a+  (»a(0),,<0>»a(0))  .  (ll) 

In  this  form  the  interpretation  of  the  energy  derivative  of  the  phase 
shift  as  the  time  delay  induced  by  the  scattering  process  is  apparent. 

I  am  indebted  to  Professor  Massey  of  University  College,  Lor- 
don,  for  his  generous  hopoitality  extended  over  part  of  the  period  during 
which  these  ideas  were  clarified. 
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A  NOTE  ON  SCATTERING  VARIATIONAL  PRINCIPLES 


by 

B.  A.  Lippmann 


ABSTRACT 


The  Schwinger  variational  principle  is  converted 
into  a  slightly  modified  form  of  the  Kohn-Hulthen 
variational  principle.  A  generalization  of  the  new 
form,  valid  for  rearrangement  collisions,  is  given. 


DIRECT  COLLISIONS 

We  start  with  the  Schwinger  variational  principle 

Tba  '  (V  V  «'a+>)  h  V  O  -  (*b_)-  I'MV|  *i+))  • 

(1) 

where 


G  = 


-  ,  IT"  ;  A  =  E  +  ie  ;  e  >  0 
o  A  -H  a  a 


(2) 


a  o 


The  last  two  terms  on  the  right  side  of  (1)  may  be  combined  to  yield: 

V  |V*a+>  +  GoV  *a+>!)  ■  -(4b'>-  <V">  |*a  +  GoV  *a+)|)  + 


(3) 


The  last  term  above  may  be  dropped  since  it  vanishes  in  the  limit  e  0. 


*  Submitted  to  ARP.,  for  publication  clearance. 
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Treating  the  first  term  on  the  right  side  of  (1)  similarly, 

(V  V  <L))  ■  (V'W  k  +  GCV  *a+>l)  -  Ic(V*a>  t4) 

where,  as  before,  the  last  term  vanishes  in  the  limit  e  -*•  0. 

As  a  consequence,  (1)  takes  on  the  new  form 


T.  =  (-f,,(A  -H  ) 
da  \  b  a  o 


$  +  G  V 

a  o  a  I  / 


(A  -H)  f  +  G  V  (+) 
\  b  a  a 


o  a 


i) 


(5) 


where  the  limit  e  -*•  0  is  understood. 

It  is  clear  that  t|/+^  and  the  quantity  in  the  { - )  are  equally 

a 

suitable  as  trial  functions  and  are  equal  for  the  exact  state  vector. 
Hence,  if  v;e  regard  { - }  as  the  trial  function,  and  symbolize  it  by 

a 

when  we  pass  to  the  limit  e  -*■  0,  a  slightly  modified  form  of  the  Kohn- 
Hulthen  variational  principle  results: 

Tb,  -  * n  -  k'1-  <e-»>  •  <6> 


where  E  is  the  common  energy  of  states  a  and  b. 


Starting  from  (1)  again, 

obtain 


and  combining  the  terms  involving 


(-) 
b  ’ 


we 


Tba  '  k'V  *b_>-  «a)  -  (<E-H>  *b">-  *a+>)  '  (7> 

The  variational  principle  (6)  d'ffers  slightly  from  the  usual  Kohn- 
(luithen  form;  the  first  term  on  the  right  is  more  conventionally  expressed 
as  (Tj,a^  trial'  present  form  automatically  insures  that  this  term 

emerges.  A  similar  remark  holds  for  (7). 

It  will  be  noted  from  (1)  that  the  presence  of  V  limits  all  integrals 
to  the  region  where  the  potential  differs  from  zero.  The  integrals  in  (6) 
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and  (7)  should  therefore  be  restricted  to  the  same  regions.  If  the  trial 
functions  are  constructed  to  represent  free  particle  motion  outside  the 
region  where  the  potential  is  different  from  zero,  this  restriction  of  the 
region  of  integration  is  automatically  accomplished  in  (6)  and  (7). 

Although  the  manner  of  derivation  of  (6)  and  (7)  from  (1)  insures 
that  the  stationary  properties  of  (1)  are  maintained  in  the  forms  (6)  and 
(7),  it  is  a  simple  matter  to  verify  directly  that  (6)  and  (7)  are  station¬ 
ary  expressions  for  T^a  when  the  state  vectors  are  varied  about  their  exact 
values. 


REARRANGEMENT  COLLISIONS 

The  foregoing  may  be  extended  to  rearrangement  collisions  by  intro¬ 
ducing  a  channel  index  (a  ■  1,2,...N)  on  the  decomposition  of  the  Hamiltonian 


H 


H  + 
o,a 


(8) 


as  well  as  the  state  vectors.  When  this  is  done,  the  variational  principle 
for  the  amplitude  scattered  from  channel  a,  state  a,  to  channel  g,  state  b, 
can  be  put  in  the  two  alternative  forms: 


8;a,a 


(E-H  )  * 
o,  p 


(+) 
a  .a 


(E-H)  <p 


(+) 

a, a 


) 


'  ((E-Ho,a>  ♦.,«)-  (<E-H>  C)  '  <9) 


COMMENTS 

The  two  forms  given  in  (6)  and  (7)  treat  the  initial  and  final  states 
unsymmetrically.  Clearly,  a  symmetrical  form  of  the  variational  principle 
may  be  obtained  by  averaging  (6)  and  (7),  and  similarly  for  the  two  forms 
given  in  (9) . 

The  trial  functions  appearing  in  these  variational  principles  cor¬ 
respond  tc  diagonalizing  the  momenta  of  the  incident  and  final  states.  If 
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the  potential  is  spherically  symmetric,  forms  of  the  variational  principles 
given  above  are  easily  obtained,  in  which  the  angular  momentum  is  diagonal, 
by  introducing  the  conventional  resolution  of  plane  waves  to  angular  mo¬ 
mentum  eigenstates.  In  the  new  forms,  the  T-raatrix  for  angular  momentum, 

£  ,  is  given  by  -  e1  1  sin  6^;  only  the  radial  parts  of  the  Hamiltonian 
operators  H  and  H  appear,  and  the  wave  functions  are  purely  radial.  To 
verify  that  the  expressions  for  each  angular  momentum  state,  £,  are  indi¬ 
vidually  stationary,  we  assume  that  the  state  vector  is  given  exactly  in 
the  plane  wave  form  for  all  £'  t  £;  varying  the  quantities  referring  to  £ 
only,  the  stationary  property  of  the  original  (plane  wave)  expression  then 
insures  that  the  terms  involving  £  alone  constitute  a  stationary  expression. 
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INTEGRAL  EQUATIONS  FOR  MULTI-CHANNEL  COLLISIONS 

by 

B.  A.  Lippman 


ABSTRACT 


Two  new  integral  equations  to  describe  multi-channel 
collisions  are  given:  one  uses  an  explicit  cut-off 
radius;  the  other  does  not.  The  new  integral  equations 
possess  symmetrical  kernels,  thereby  facilitating  the 
application  of  the  Pauli  principle  "and  permitting  the 
asymptotic  scattered  amplitude  in  each  channel  to  be 
derived  simply.  The  limit  e  -*•  0  does  not  result  in 
a  "homogeneous"  integral  equation. 


The  integral  equations  appropriate  to  each  of  the  channels  in  a 
multi-channel  collision  are  here  unified  by  the  use  of  channel  projection 
operators.*  Two  new  integral  equations  result  corresponding  to  the  choice 
of  two  different  types  of  channel  projection  operators  used:  one  type  uses 
an  explicit  cut-off  radius  and  the  other  does  not. 


The  new  integral  equations  treat  all  channels  on  an  equal  basis; 
permit  simple  derivations  of  the  asymptotic  form  of  the  scattered  amplitude 
in  each  channel  —  from  which  the  transition  matrix  elements  -'.re  easily 
identified;  possess  symmetrical  kernels,  thereby  facilitating  the  applica¬ 
tion  of  the  Pauli  principle;  and  eliminate  the  difficulties  in  principle  aris- 

2 

ing  from  the  presence  of  "homogeneous"  equations  in  some  of  the  channels. 

*  Submitted  to  ARPA  for  publication  clearance. 

1  Related  operators  have  previously  been  introduced  by  H.  Feshback, 

Ann.  Phys.  (N.Y.)  Jj,  357  (1958);  1J9,  287  (1962);  and  B.  A.  Lippmann, 
Lockheed  Technical  Report  NONR  3368(00) ,  1964  (unpublished) . 

2  This  apparent  lack  of  uniqueness  in  multi-channel  collision  theory 
has  led  to  the  introduction  of  the  Fadeeyev  equations.  There  is  no 
lack  of  uniqueness  if  e  ^  0. 
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The  derivation  is  Riven  in  a  general  form  and  also  illustrated  for 
the  case  of  electron  scattering  by  atomic  hydrogen.  The  electrons  are 
first  considered  distinguishable:  the  Pauli  principle  is  taken  into  account 
later. 


For  electron-H  scattering,  there  are  four  channels.  Assuming  shield¬ 
ed  Coulomb  potentials,  channel  1  corresponds  to  electron  1  moving  freely 
and  electron  2  interacting  with  the  proton;  channel  2  corresponds  to  elec¬ 
tron  2  moving  freely  and  electron  1  interacting  with  the  proton;  channel  3 
corresponds  to  both  electrons  moving  freely;  channel  4  corresponds  to  both 

electrons  interacting  with  the  proton  and,  if  convenient,  also  with  each 
3 

other.  There  is  some  arbitrariness  in  the  division  of  the  interactions 

between  H  and  V  ,  where  H  =  H  +  V  (a  =  1,2, 3, 4). 
o,a  a  o,a  a 


The  basis  vectors  in  each  channel  are 


*a,l 

(1,2) 

=  ^a' 

(2) 

Xa" 

(1) 

^a.2 

(1,2) 

(1; 

xa" 

(2) 

^a,  3 

(1,2) 

=  V 

(1; 

xa" 

(2) 

^a,4 

(1,2) 

(1  j 

Ua" 

(2) 

(1) 


where  the  x’s  cor~espond  to  .  particle  wave  functions  and  the  remain¬ 
ing  symbols  describe  interacting  wave  functions.  Each  tf>  above  is  an 
eigenfui.ction  of  the  corresponding  H 

0,0 


Assuming  the  incident  wave  is  in  channel  1,  the  integral  equation 
for  the  total  wave  function  is 


+ 


G  i  vi 

o,l  1 


(+) 

a,  1 


(2) 


3  In  particular,  channel  4  can  be  defined  to  include  a  zero  width 
resonant  state,  with  being  responsible  for  the  energy  shift  and 
lifetime. 
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where 


o,l 


A  -  H  . 
a  o,l 


;  A  =  E  +  ie  ; 
a  a 


e  >  0  .  (3) 


In  channels  2,3,  and  4,  we  have  (in  the  limit  e  -+•  0) 

4<(+j  =  G  V 

a, 1  o,2  2  a, 1 

=  G  -  V,  ip(+J 
O  y  3  3  (1)  1 


G  .  V,  rp(+?  • 

o,4  4  a,l 


(4) 


The  integrations  in  coordinate  space  consist  of  the  direct  products 
of  the  individual  coordinate  spaces  of  electrons  1  and  2,  which  we  divide 
into  channels  as  follows.  In  each  space,  we  introduce  a  distance  R  around 
the  target.  Let  Rj(l)  be  the  operator  that  projects  onto  the  space  of 
particle  1  interior  to  R;  similarly,  let  Rg(l)  project  onto  the  space  of 
particle  1  exterior  to  R.  These  are  the  "internal"  and  "external"  regions 
of  1.  For  any  value  of  the  total  energy,  R  can  be  chosen  large  enough  so 
that  if  one  examines  the  wave  functions  of  the  bound  states  energetically 
possible  in  each  channel,  these  wave  functions  are  essentially  contained 
within  R. 


Vie  define  four  channel  projection  operators  corresponding  to  the  four 
channels  that  are  possible  in  this  special  case: 


C1  =  RI(2)  RE(1) 

C3  =  RE(1)  V2) 


C2  =  RI(1)  R£(2) 

C4  =  Rj(l)  Rj (2) 


(5) 


The  channel  operators  are  mutually  exclusive  and  their  sum  covers  the  direct 
product  space  completely; 


cicj  ■  ci  s«  =  2  ci  •  1 

i 


4  B.  A.  Lippmann,  Phys.  Rev.  102,  264  (1956). 


(6) 
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We  now  multiply  Eq.  (2)  by  C^,  and  the  successive  lines  of  (A)  by 
C^,  C 2*  and  C^,  and  add,  thereby  obtaining  the  first  of  the  integral  equa¬ 
tions  we  have  been  seeking 


C1  *a 


C  G  V 
a  o ,  u  u 


a 


(7) 


In  the  limit  t  ■*  0,  the  second  term  on  the  right  yields  the  scattered 
amplitude  in  each  channel.  That  is,  the  projection  operators  effect  a 
separation  between  those  states  which  only  exist  in  the  region  internal  fo 
R  and  those  which  exist  in  the  external  region  as  welx:  these  are  the  dis¬ 
crete  and  continuum  states  respectively,  if  R  has  been  chosen  large  enough, 
the  discrete  states  in  the  internal  region  that  are  energetically  accessible 
maintain  their  full  c’-thogonality  properties  unimpaired.  The  continuum  states 
in  the  external  region  that  combine  with  these  bound  states  to  conserve 
energy  are  unaffected  by  starting  the  integration  at  a  finite  value  of  R; 
the  continuum  states  that  do  not  conserve  energy  are  exponentially  damped 
in  the  external  region.  Thus,  as  R  approaches  infinity,  these  exponentially 
damped  states  contribute  nothing.  For  the  "propagating"  states  in  the  ex¬ 
ternal  region,  R  may  be  permitted  to  approach  zero. 


introducing  the  basis  vectors  for  each  channel  explicitly,  (7)  be¬ 
comes 


C1  *a,l+  2  I 

a  b 


C 

a 

A 

a 


T 


b,ct;a,l 


(8) 


where  the  transition  amplitudes  are  given  by 


D  |  3  )  Cl 


*(+> 
cl  9  J. 


(9) 


We  have  assumed  distinguishable  particles  in  deriving  (8).  However, 
if  the  particles  are  indistinguishable,  (7)  may  be  adapted  to  reflect  this 
explicitly  by  multiplying  both  sides  by  the  proper  symmetr izat ion  or  anti- 
symmetrizat.ion  operator,  s.  With  the  limit  l  ►  0  understood,  (8)  is  then 
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replaced  by 


3  C1  *3,1  +  S  S 

a  b 


C 

a  b,a 

A  -  E, 
a  b 


V 

a 


3  III 


(+) 

£,1 


The  sum  in  (10)  includes  states  that  differ  only  by  an  interchange  of 
particles  1  and  2.  If  the  summation  is  restricted  so  that  these  states  are 
excluded  —  that  is,  the  summation  is  only  over  "  'istinct"  channels  —  (10) 
becomes 


■  <+> 
a,l 


s  C.  * 


1  a,  1 


I 


s  C  4> 


a  b,a 


a,b 

"distinct" 

states 


A  - 

a 


b.a 


V  s 
a 


(+)  1 
'a.l/ 


(ID 


Thus  the  transition  amplitude,  when  the  Pauli  principle  is  operative, 
is  given  by 


T. 


B;a,a 


V  * 


(12) 


The  Drojection  operators  defined  so  far  correspond  to  a  definite 
choice  of  a  cut-off  radius  However,  (8)  also  holds  for  the  channel  pro¬ 
jection  operators,  defined  below,  which  are  independent  of  the  cut-off 
radius. 


These  Ca  are  defined  in  the  terms  of  the  particle  energies  as  follows. 
In  each  channel,  say  "a,"  the  corresponding  H  defines  the  enetgies  of 

O  f  Cl 

particles  1  and  2  —  when  one  of  them  is  in  the  asymptotic  region  —  unam¬ 
biguously.  We  can  therefore  redefine  the  C  : 


n(-E2)  n(E1) 


n(-E1)  n(E2) 


n(E^)  n(E2) 


=  n(-E1)  n(-E2) 


2r,(x) 


1  + 


R  * 


tl3) 
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The  previous  considerations  hold  unchanged  for  the  new  definitions 
of  the  channel  projection  operators,  in  particular,  (8)  no  longer  requires 
the  introduction  of  a  cut-off  radius. 
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CALCULATION  OF  DISSOCIATIVE 
ATTACHMENT  IN  HOT  02* 

by 

T.  F.  O'Malley 


ABSTRACT 


The  very  striking  temperature  dependence,  of  the  dis¬ 
sociative  attachment  cross  sections  found  experimentally 
in  O2  has  been  reproduced  by  a  semi-empirical  calcula¬ 
tion,  assuming  a  Maxwell  distribution  of  vibrational 
(v)  and  rotational  (r)  states.  A  recently  derived  ex¬ 
pression  for  the  cross  section  ov^r  is  used,  and  the 
final  state  potential  curve  of  Of  is  parametrized  to 
fit  the  experiments.  The  temperature  shift  is  found  to 
be  caused  by  the  effect  on  excited  vibrational  states 
of  the  rapidly  varying  "survival  probability,"  a 
measure  of  the  competition  between  autoionization  and 
dissociation.  By-products  of  the  calculation  are  the 
potential  curve  for  the  final  dissociating  state  of  O2- 
together  with  an  estimate  of  its  autoionization  width, 
ra(R),  and  also  the  contribution  of  this  state  to  the 
total  cross  section  for  electron  energy  loss  to  vibra¬ 
tional  excitation  and  dissociation. 


It  has  been  found  experimentally  by  Fite  et  al.  that  the  cross 
section  for  dissociative  attachment  (DA)  of  electrons  to  02  has  an  extra¬ 
ordinarily  large  dt  •mdence  on  the  temperature  of  the  C^.  The  shift  and 
broadening  with  temperature  waj  such  that  at  .  '',)°K  the  cross  section  peak 
was  shifted  to  lower  energy  by  1  eV  while  the  apparent  onset  v.\-  s  reduced 


*  To  appear  in  Pnys.  Rev.,  March  1967. 

^  W.  L.  Fite  and  R.  T.  Brackmann,  Proc.  4th  International  Conference 
on  the  Ionization  Phenomena  in  Gase3,  Paris,  1963,  Vol  I,  p.  21. 

W.  L.  Fite,  R.  T.  Brackmann,  and  W.  R.  Henderson,  IVth  International 
Conference  on  the  Physics  of  Electronic  and  Atomic  Collisions,  Science 
Bookcraf ters,  Inc.,  Hastings,  N.Y .  (1965),  p.  100. 
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by  over  2  eV.  The  present  calculation  was  done  in  an  attempt  to  provide 
some  understanding  of  this  remarkable  effect. 


It  is  first  assumed  that  the  direct  effect  of  the  temperature  (T)  of 

the  0o  is  to  produce  a  Maxwell  distribution  of  vibrational  (v)  and  rota- 

t'onal  (r)  states.  The  cross  section  o(T,£),  where  E  is  the  electron 

energy,  is  then  the  Boltzmann  average  of  the  cross  sections  a  (£)  from 

v»  r 

each  of  the  individual  excited  states,  thus 


(T,E) 


v=v  . 


mm 


s 


-(E  +E  )/kT 
v  v  r 
Ne 


o 


v,r 


(E) 


r=r  . 


min 


(1) 


where  N  is  the  boltzmann  no-malizat ion  factor,  and  E  and  E  are  the  vibra- 

v  r 

tional  energy  and  rotational  energies  of  the  0^  molecules,  respectively. 

The  quantities  v  .  and  r  .  are  subject  to  the  threshold  requirement  for 
min  mm  J 

the  process,  viz.  E  +  E  +  E  >  3.6  eV,  and  the  cross  section  o  is 

v  r  —  v,r 

given  by  the  theoretical  expression^ 


a 

v,  r 


(E) 


4n 


r 


-P 


(2) 


The  F  ,  F  and  p  as  wel1  as  x  are  evaluated  at  the  final  state  turning 

<i  Cl  V 

point  Rj.(E,v,r).  The  exact  definition  of  the  quantities  is  given  in  Ref.  2. 

In  particular,  k  is  the  electron's  momentum,  g  is  an  angular  momentum  and 

spin  factor.  F  and  T  are  the  total  and  partial  (for  the  ground  state 
a  a ,  x 

X)  aatoionization  widths,  '.e.,  T  =  Z.  F  .  F.  is  a  width  for  dis- 

a  j  a,j  d 

sociation,  x  is  the  vibrational  wave  function,  R  the  turning  point  is  in 

V  b 

vibrational  units  and  p  is  twice  tiie  imaginary  part  of  the  final  state 

_  -V  <1 

phase  shift  f  .  The  factor  e  ^  is  called  the  survival  probability  or 
survival  factor.  It  is  given  approximately  by 


exp 


R 


c 

r 

-  I  F  (R)dR/hv(R) 

J  a 

=  exp 

-  j  ‘  a(R)  dt(R)  /  1i 

■ 

(3) 


T.  F.  O'Malley,  Phys.  Rev.  JJ>0,  14  (1966). 
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where  R  is  the  point  where  i'  goes  to  zero,  v  is  the  velocity  of  dissocia- 

C  3 

tion  and  t  the  corresponding  classical  time.  This  factor  results  from  the 

competition  between  autoionization  and  dissociation.  Where  the  above  ex- 

3 

pression  is  not  sufficiently  accurate,  an  exact  expression  can  be  used. 

The  effect  of  rotation  in  (2)  is  implicit.  It  adds  a  centrifugal  term  to 
the  potential  curves  for  initial  and  final  vibrational  motion,  and  at  low 
energies  it  helps  put  the  total  energy  above  the  threshold  and  thus  de¬ 
termines  v  .  in  Eq.  (1). 
rain 

The  cross  section  (2)  is  completely  determined  when  the  potential 

energy  curves  V,  and  V,  are  specified  including  the  width  T  of  the  latter, 
it  a 

which  is  a  resonance  or  autoionizing  state.  Since  V  ,  the  ground  state 

4  1 

curve  for  0^  is  well  known,  it  is  necessary  only  to  know  the  final  curve 
together  with  its  width  T  in  order  to  determine  the  entire  family  of 
cross  sections  given  by  (2).  Accordingly,  the  procedure  adopted  in  the 
present  work  was  to  parametrize  this  potential  curve  and  its  width.  The 
parameters  are  then  chosen  with  a  view  to  fitting  the  experimentally  de¬ 
termined  a (T, E) ,  in  particular  as  given  by  Ref.  1,  to  the  extent  that  this 

is  possible.  The  curve  V  is  first  expanded  in  a  series  about  R  =1.21 

t  o 

(the  0^  equilibrium  distance) ,  as 

vf(R)  =  Eq  -  V'AR  +  V"AR2  -  (1/6)  v"  *AR3  (4) 

where  AR  =  R  -  R  .  For  this  width,  we  mav  either  parametrize  T  (R)  and 
o  'a 

then  derive  p  frpm  Eq.  (3)  or  a  more  exact  formula,  or  else  parametrize 
p(R)  directly  and  derive  T  .  The  latter  course  proved  easiest,  ar was 
done  in  the  form 

p(R)  =  p'  (Rc  -  F)n  (R  £  Rc)  (5) 

3  J.  N.  Bardsley,  A.  Herzenberg,  and  F.  Mandl,  Atomic  Collision  Pro¬ 
cesses  (M.R.C.  McDowell,  ed.),  North  Holland,  Amsterdam,  1964,  p.  415. 

^  F.  R.  Filmore,  J.  Quant.  5pectry.  Rad.  Trans.  _5,  369  (1965). 
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and  p  =  0  beyond  Rc«  Lest  the  form  of  (5)  appear  arbitrary,  this  is 
actually  the  exact  form  which  follows  from  (3)  in  the  limit  that  the  po¬ 
tential  curve  V  is  approximately  linear  between  R_  and  R  and  T  has  its 
1  L+l/2  k  c  a 

threshold  behavior  P^  “  (R^  -  R)  ,  where  L  is  the  electron  s  asymptotic 

angular  momentum.  In  this  case  n  would  be  equal  to  L  +  1.  However,  the 

present  results  cover  a  wider  range  and  so  (5)  represents  simply  a  3-pa- 

rameter  approximation  to  the  true  behavior.  R£  is  interpreted  as  the 

crossing  point  between  V  and  the  target  potential  curve  into  which  it 

autoionizes,  which  is  where  the  partial  autoionization  width  P  .  and 

a*  J 

therefore  its  contribution  to  p  goes  to  zero  in  the  Born-Oppenheimer  ap¬ 
proximation  . 

It  was  found  that  the  best  fit  to  the  experimental  results^  was 
obtained  with  the  parameters  in  (4)  and  (3)  chosen  as  follows: 

E  =  7.8  i  0.2  eV,  V'  =  27  i  2  eV/X,  V"  -  110  t  20  eV/A2  , 
o 

V"  '  =  300  i  100  eV/A3  (6) 

and 

o- i / ?  o 

n  =  1.5  £  0.5,  p'  =  41  t  3  A  ,  R  =  1.44  i  0.03  A  .  (7) 

c 

The  uncertainties  expressed  for  the  various  quantities  ate  meant  to  indi¬ 
cate  the  sensitivity  of  the  fit  to  the  various  parameters,  and  do  not 
necessarily  reflect  any  judgment  as  to  absolute  errors. 

The  resulting  cross  sections  at  300  and  at  2100°K  are  plotted  in 
Fig.  1  together  with  the  experimental  results  at  these  temperatures.  The 
theoretical  curve  for  3000°  is  also  shown.  The  normalization  of  the  300 
cross  section  was  taken  mainly  from  the  work  of  Schulz.3  The  theoretical 
curves  in  Fig.  1  have  been  averaged  over  a  broad  distribution  of  electron 
energies  chosen  to  approximate  that  employed  in  the  experiment 

o  (T,E)  =  /  dE'  f(E')  n(T,E' )  (8) 

av 

^  G.  J.  Schulz,  Phys.  Rev.  128,  178  (1962).  See  also  D.  Rapp  and 
D.  D.  Briglia,  J.  Chem.  Phys.  4^3,  1480  (19b5). 
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where  the  electron  distribution  function  f(E')  was  taken  to  be  a  shifted 
Maxwellian  corresponding  to  about  2800  K.  The  experimental  2100  cross 
section  represents  the  experimental  signal  normalized  to  agree  in  magnitude 
with  the  presently  calculated  2100  cross  section,  which  is  1.3  times  as 
large  as  that  at  room  temperature.  Though  the  magnitude  of  the  cross 
section  (as  opposed  to  the  signal)  is  not  given  explicitly  in  Ref.  1,  the 
ratio  of  1.3  found  here  is  consistent  with  the  cross  sections  implied  by 
the  discussion  in  these  papers. 

The  theoretical  fit  to  the  experiments  may  be  seen  from  Fig.  1  to  be 
good,  except  that  it  falls  somewhat  belcw  at  the  lowest  electron  energies. 
The  experimental  plateau  above  8  eV  is  apparently  due  to  some  other  cause. 
In  Table  1,  the  calculated  cross  sections  o(T,E)  are  tabulated  for  a  num¬ 
ber  of  temperatures.  Note  that  unlike  Fig.  1,  these  are  not  averaged  over 
the  electron  distribution  of  Eq.  (8). 

In  the  course  of  the  fitting  process,  it  was  found  that  the  three 
parameters  Eq,  V'  and  V"  of  the  potential  curve  are  essentially  de¬ 
termined  by  the  room  temperature  data,  except  that  E  is  raised  by  the 

—  p 

presence  of  e  .  The  quantity  V’"  does  not  affect  the  300°  result  and  was 
chosen  mainly  for  its  influence  on  D(E)  (see  Table  2).  The  expansion  (4) 
of  Vj  should  hold  till  roughly  R  =  1.55  $. 

—  * 

The  final  state  potential  curve  for  the  system  0^  is  shown  in 

Fig.  2,  together  with  other  relevant  states  of  09  and  0?  taken  from 

A  ^  ^ 

Gilmore's  curves.  It  is  believed  that  this  potential  curve  together  with 

the  approximation  (5)  and  (7)  to  the  width  (see  also  footnote  6)  consti¬ 
tute  a  significant  by-product  of  the  present  work.  Note  that  the  curve 


0  Strictly  speaking,  one  should  put  at  least  two  terms  like  (5)  into 
one  for  each  set  of  states.  This  was  tried  and  found  to  give 
essentially  the  same  cross  sections.  The  parameters  were  found  very 
roughly  to  be  (with  n  =  1.5)  R^  =  1.52,  r'j  -  lb;  Rp7  =  1.39, 
r ' 7  =  27.  This  gives  approximately  Tj  =  1.5  (1.52-R)  and  1^  =  2.6 
(1T39-R),  while  Tx  £  0.6  r ^ .  However,  the  relative  magnitude  of 
and  .2  was  found  to  be  somewhat  arbitrary  without  additional 
experimental  information. 
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TABLE  1 

Dissociative  attachment  cross  sections  a(T,E)  from  Eq.  (1) 

-21  2 

at  various  temperatures,  in  units  of  10  cm  .  The 
electron  energy  E  is  in  eV  and  T  in  °K.  A  sharply 
defined  electron  energy  is  assumed. 


E/T 

2600° 

2200° 

1800° 

1400° 

300° 

1.0 

5 

1 

0 

0 

0 

1.5 

22 

4 

0 

0 

0 

2.0 

98 

27 

4 

0 

0 

2. 5 

312 

119 

29 

3 

0 

3.0 

650 

306 

104 

20 

0 

3.5 

1070 

610 

271 

77 

0 

4.0 

1^30 

1010 

560 

224 

1 

4.5 

1?20 

1430 

940 

500 

21 

5.0 

2120 

1750 

1340 

890 

155 

5.5 

2070 

1870 

1610 

1280 

520 

6.0 

1810 

1750 

1650 

1500 

1040 

6.5 

1420 

1440 

1450 

1440 

1380 

7.0 

1000 

1G50 

1100 

1160 

1310 

7.5 

630 

670 

720 

770 

930 

8.0 

360 

381 

405 

434 

520 

8.5 

183 

190 

197 

206 

229 

9.0 

82 

■ 

83 

83 

83 

83 
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TABLE  2 

lable  of  activation  energies,  D(E),  experimental  and  theoretical 

in  eV,  given  in  different  approximation.  D's  are  in  eV. 

D  comes  from  o  of  Eq .  (8),  i.e.,  with  electron  spread, 
av  av 

D  corresponds  to  a  with  no  electron  spread. 

C  is  the  corresponding  coefficient  from  Eq.  (10) 
in  units  of  10  cm  . 

was  compared  with  rotation  entirely  neglected  (vibration  only). 


E 

D 

exp 

D 

av 

D 

D 

V 

C  (E) 

1.0 

— 

1.59 

2.43 

2.61 

370 

1.5 

— 

1.36 

2.05 

2.12 

280 

2.0 

1.0 

1.14 

1.64 

1.66 

200 

2.5 

.92 

.93 

1.18 

1.18 

72 

3.0 

.79 

.72 

.89 

.91 

39 

3.5 

.58 

.55 

.67 

.68 

23 

3.75 

— 

.47 

.59 

.58 

19 

4.00 

.29 

.40 

.50 

.50 

15 

4.25 

— 

.34 

.42 

.42 

12 

4.50 

.15 

.28 

.35 

.35 

10 

4.75 

— 

.22 

.28 

.28 

8 
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goes  below  the  dissociation  limit  at  about  R  =  1.44  A  and  becomes  attrac- 

tive.  T'  j  dotted  line  from  about  R  =  1.35  to  2.1  %  represents  the  way  in 

which  this  curve  is  assumed  to  connect  up  with  the  known  asymptotic  form 
2  ,  -* 

of  the  ti  curve  for  0^  •  Although  no  use  is  made  of  this  dotted  region 

of  the  curve,  the  configuration  is  fully  consistent  with  the  assumed 
electronic  structure  of  the  state ,  which  would  consist  of  an  electron  attach¬ 
ed  loosely  to  a  combination  A  and  C  core.  It  follows  that  this  state  probably 
has  a  dissocial  ion  energy  of  the  order  of  1  -j  eV. 

The  heuristics  of  choosing  the  parameters  for  p  in  (7)  may  be  of 
interest.  They  were  chosen  to  fit  the  2100°  cross  section.  The  primary 
feature,  namely  the  shift  of  the  peak  by  1  eV,  was  easily  reproduced  by 
simply  varying  the  amplitude  p'.  However,  the  obvious  choice  of  n  =  1/2 
(corresponding  to  a  constant  F  )  gave  a  cross  section  that  was  either 

3 

nearly  square  in  shape  or  double-peaked,  depending  on  R^.  It  was  only 
when  n  was  varied  to  between  1  and  2  that  bell-shaped  curves  resulted. 

The  value  of  1.5  which  fits  best  is  gratifyingly  close  to  the  assumed  cor¬ 
rect  threshold  value,  which  as  mentioned  would  give  n  =  2.  Finally,  fitting 
the  observed  cross  section  width  fixed  R^  loosely  at  1.44  i  0.03  X.  The 
value  1.44  for  Rc  was  a  little  surprising.  Reference  to  Fig.  2  shows  that 
this  lies  squarely  between  the  crosring  of  the  A,C  system  at  R  •'  1.39  and 
the  crossing  of  the  x»a,b  states  at  R  >_  1.55.  R^  had  been  expected  to 
coincide  with  one  or  the  other  of  these  crossings.  The  intermediate  value, 
somewhat  closer  to  A  and  C,  is  Interpreted  as  meaning  that  there  is  con¬ 
siderable  autoionization  to  all  of  the  5  allowed  states,  with  the  largest 
part  going  to  A  and  C.  The  value  of  1.44  is  therefore  taken  as  an  average 
or  compromise  value. 

The  foregoing  results,  (6)  and  (7),  were  derived  with  the  cross  sec 

JO 

tions  arbitrarily  normalized.  Taking  the  300  normalization  of  1.4  x  10 
cm~  from  Ref.  5  allows  F  in  (2)  to  be  determined.  First,  assuming  the 

i  a  * x 

0^  state  to  be  i!^,  it  follows  that  the  weighting  factor  g  equals  2/3. 
Further,  it  follows  from  n  =  1.5  in  (7)  that  the  width  F  is  roughly  linear 

3 

in  electron  energv.  Substituting  all  this  into  Eq .  (2)  with  v  =  0  and  c 
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small  yields  for  the  partial  capture  width  from  the  stat.  X 

P  „  =  0.034  E  . 

X 


(9) 


This  was  Incorporated  In  Fig.  1  and  Table  1.  The  total  width,  T  ,  also 

can  be  computed^  from  p  as  determined  by  (5)  and  (7).  As  T  is  not  used 

directly,  this  was  done  only  roughly,  in  a  linear  approximation,  with  the 

result  F  (R)  %  4  (1.44  -  R) .  Near  the  cross  section  peak  then  T  is 
a  a,  a 

somewhat  more  than  1/4  of  I'  ,  which  Is  also  believable. 

a 


An  analysis  of  the  physical  content  of  the  ca’  alation  is  made  in 
Fig.  3,  where  the  2100°  cross  section  is  plotted  in  five  successive  approxi¬ 
mations.  In  curve  a,  we  set  v  =  r  =  p=  0.  This  would  correspond  to  the 
low  temperature  cross  section  without  p.  Including  excited  vibrational 
states  in  (b)  broadens  the  curve  considerably.  Applying  the  rapidly  vary¬ 
ing  survival  factor  e  P  shifts  (b)  to  (c) ,  while  the  narrower  curve  a  is 
shifted  only  to  c1.  (Differences  in  magnitude  are  suppressed  in  Fig.  3.) 

To  a  good  approximation  for  O2,  the  difference  between  c  and  c'  represents 
the  entire  temperature  effect.^  If  p(E)  or  its  energy  variation  were 
small,  as  for  some  other  molecules,  then  the  curves  a  and  b  would  describe 
the  effect  and  there  would  be  broadening  but  no  shift.  Curves  d  and  e  , 
representing  0  and  a  respectively,  show  the  modest  shift  caused  by  rota¬ 
tional  states  (agreeing  with  the  analysis  of  Ref.  2)  and  the  shift  and 
tailing  effect  of  a  broad  electron  distribution.  It  can  be  seen  from  (d) 
that  excited  rotational  states  are  of  negligible  importance  compared  to 
vibrational  states. 


Figure  4  analyzes  the  calculation  from  a  slightly  different  point  of 

view.  The  cross  sections  ov  for  attachment  frem  the  vth  state  are  plotted 

individually.  (The  small  dependence  on  r  is  here  suppressed.)  For  each 

increase  in  v,  there  is  a  large  shift  to  the  left  and  a  subs  antial  increase 

in  magnitude  (o_  is  20  times  as  large  3s  o  ) .  Note  that  each  cross  section 
J  o 

is  cut  off  below  its  threshold  (E  +  E^  =  3.6  eV) .  These  curves  show 
graphically  the  basis  for  the  observed  shift  and  the  large  low  energy  tail. 

7  This  is  in  agreement  with  the  qualitative  suggestion  of  Y.  N.  Demkov, 

Phys.  Letters  15.,  235  (1965).  See  also  Ref.  2. 
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A  very  interesting  aspect  of  the  work  of  Fite  et  al.  was  the  finding 
that  at  the  lower  energies  (E  5  eV  and  T  \  1000  K)  the  cross  sections 
follow  very  closely  the  formula 


o(T,E)  =  C(E) 


-D(E)/kT 


(10) 


Exactly  toe  same  was  found  for  the  presently  computed  results.  In  Fig.  5, 
log  o  is  plotted  against  1/T  and  the  data  are  seen  to  lie  on  straight  lines 
from  roughly  800  to  3000'K.  The  slopes,  D(E),  which  were  called  activation 
energies,  have  a  clear  interpretation  *n  the  present  work.  They  are  the 
average  internal  energy  (vibrational  plus  rotational)  of  the  molecule  which 
maximizes  the  summand  in  Eq.  (1) ,  showing  the  internal  states  which  contribute 
most  strongly  to  the  temperature-averaged  cross  section.  In  columns  1  and  2 
of  Table  2,  the  experimental  values  Dexp(E)  are  compared  with  the  theoreti¬ 
cal  values  D  (E)  taken  from  c  of  Eq.  (8).  The  agreement  is  seen  to  be 
av  av 

within  0.1  eV  everywhere.  (Probably  the  discrepancy  at  the  largest  E’s 
would  be  reduced  if  the  experimenters  had  plotted  a  rather  than  the  uncor¬ 
rected  current.)  The  third  column  shows  D(E) ,  the  present  result  with 
sharply  defined  electron  energies  from  Eq.  (1).  The  very  large  difference 

between  D  and  D  is  due  to  the  tailing  effect  of  the  broad  electro. i  dis- 
av 

tribution.  In  the  fifth  column,  the  corresponding  coefficient  C(E)  from 
Eq .  (6)  is  listed.  The  fourth  column,  labeled  Dv>  was  calculated  with 
vibrational  states  only  (no  rotational  excitation) .  This  corresponds  to 
the  data  plotted  in  Fig.  6,  which  differs  very  little  from  Fig.  5.  The 
difference  between  D  and  is  some  measure  of  the  relative  importance  of 
rotation,  which  seems  to  be  very  small  except  below  about  2  eV.  This  is 

In  disagreement  with  the  assessment  of  Fite  et  al. ^  In  this  connection, 

o 

it  was  found  that  at  2100  rotational  states  with  E  >>  kT  never  contributed 

r 

significantly  to  the  sum  in  Eq.  (1)  (unlike  the  situation  with  the  vibra¬ 
tional  states),  until  E  fell  below  about  2  eV.  Ever  at  the  lowest  energy 
of  1.0  eV  the  average  excitation  of  the  rotational  states  contributing  to 
(1)  was  not  much  more  than  20%  of  that  of  the  average  vibrational  state 
(v  13). 
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Figure  5.  Log  o(T,E)  fr™"  the  present  calculation  is  plotted  vs  1/T  for  va 
electron  energies  of  5  eV  and  below.  The  empirical  formula  (10) 
resents  straight  lines  on  this  plot.  Like  the  experimental  resu 
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Figure  6.  The  same  as  Fig.  5  but  from  a  calculation  do  le  without  excited  rota¬ 
tional  states.  The  straight  line  behavior  indicates  that  the  relation 
given  by  Eq.  (10)  does  not  necessarily  depend  on  rotational  excitation. 


UNCLASSIFIED 


It  is  interesting  in  Table  2  that  for  electron  energies  well  below 
the  threshold  of  3.6  eV,  the  values  o'  D(E)  (not  D  )  converge  more  or 
le3S  to  3.6-E,  the  internal  energy  necessary  to  raise  the  system  above  the 

o 

threshold.  This  would  indicate  that  the  results  for  very  high  vibrational 
states  are  no  longer  sensitive  to  the  details  of  the  wave  functions,  but 
that  these  may  be  replaced  by  a  suitable  average  value,  in  the  statistical 
spirit.  Thus,  the  simple  concept  of  an  activation  energy  does  become  ap¬ 
propriate  in  this  limit.  It  seems,  however,  that  at  the  very  lowest  energies 
(below  2  eV)  the  dynamical  effects  of  rotation  tend  to  bring  D(E)  (as 
opposed  to  Dv)  below  the  simple  threshold  value  of  3.6-E. 

Finally,  the  consequences  of  the  present  analysis  for  processes  other 

than  dissociative  attachment  should  be  pointed  out.  Consider  only  the  low 

3 

temperature  cross  section  for  simplicity.  The  accepted  interpretation  of 
Eq.  (2)  is  that  the  DA  cross  section  is  equal  to  a  capture  cross  section  to 

-*  -p 

the  state  of  02  times  e  ,  the  chance  of  surviving  to  dissociate  befc  e 
autoionizing  to  some  vibrational  level  of  one  of  the  I  -e  lowest  0^  states. 
Given  the  presently  determined  parameters  (6)  and  (7),  it  follows  that  the 
capture  cross  section  e+P  a  ,  is  roughly  70  times  that  for  DA,  rith  a 

L»A  2 

peak  value  of  about  1.0  x  10“ 16  cm  at  7.8  eV.  From  R  *  1.44  we  would 

9  c 

expect  that  at  least  half  of  this  capture  cross  section  will  decay  to  the 
states  A  and/or  C,  resulting  in  dissociation,  while  the  remainder  auto- 
ionizes  to  the  states  X,  a  and  b,  (about  half  of  it  to  X),  resulting  in 
excitation  of  most  of  the  vibrational  levels  of  these  states.  Whether  the 
contribution  of  this  to  the  respective  cross  sections  is  significant  is 
not  clear  at  present. 


8  This  relation  was  obtained  previously  by  M.  A.  Fineman  from  similar 
experiments  done  by  J.  W.  McGowan.  (Private  communication.) 

^  See  also  Footnote  6. 
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QUASl-ADIABATIC  MOLECULAR  STATES:  ^ 

AN  ALTERNATIVE  TO  THE  NON-CROSSING  RULE 

by 

Thomar.  F.  O'Malley 


ABSTRACT 

Because  of  the  non-crossing  rule,  the  conventional 
definition  of  adiabatic  states  has  been  increasingly 
regarded  as  unsuitable  for  describing  collisions. 

An  alternative  definition  based  on  a  generalization 
of  the  resonance  is  given,  defining  states  which  may 
cross  freely.  As  an  immediate  application,  these 
states  tend  to  resolve  a  serious  discrepancy  in  the 
interpretation  of  recent  He*  recombination  data. 

The  non-crossing  ruin  (NCR)  of  von  Neuman  and  Wigrier1  which  states 

that  molecular  states  of  the  same  symmetry  cannot  cross,  has  been  a  source 

2 

of  growing  dissatisfaction  to  many  molecular  physicists  interested  in 
atomic  collisions  and  reactions.  Nevertheles  **  he1'  been  accepted,  al¬ 
though  sometimes  -  3luctantly,  as  the  apparent  price  that  one  must  pay  for 
having  an  "adiabatic"  definition  of  electronic  states,  i.e.,  a  definition 
in  which  the  Internuclear  coordinate  R  is  held  temporarily  fixed.  It  is 
the  purpose  of  the  pres  mt  note  to  point  out  that  there  exists  an  accept- 


*  Submitted  to  Phys.  Rev.  letters. 

^  J.  Von  Neumann  and  E.  P.  Wigner,  Phys.  Zeits.  3£,  467  (1929). 

^  See  for  example  W.  L.  Lichten,  Phys.  Rev.  131 ,  229  (1963).  Lichten's 
emphasis  on  the  physical  importance  of  certain  pseudocrossing  single 
configuration  states  is  very  much  in  the  spirit  of  the  present  work. 
These  states  would  be  a  first  rough  approximation  to  the  4>r  of  the 
present  paper.  F.  T.  Smith,  D.  C.  Lorents,  W.  Aberth  and  R.  P.  Marchi, 
Phys.  Rev.  Letters  _15_,  >42  (1965)  have  been  forced  to  assume  similar 
crossing  curves.  See  also  W.  R.  Thorson  and  S.  A.  Boorstein,  in  IVth 
International  Conference  on  the  Physics  of  Electron  and  Atomic  Col¬ 
lisions  (Science  Bookcraf ters ,  Hastings,  N.Y. ,  1964),  p.  218. 
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able  and  precise  definition  of  adiabatic  electronic  states  which  are  not 
constrained  by  an  NCR  and  which,  as  a  result,  are  able  to  approximate  the 
physical  states  which  are  actually  followed  in  various  types  of  rearrange¬ 
ment  collisions. 

The  simplest  way  of  presenting  the  NCR  is  in  terms  of  two  approximate 
trial  functions,  4^  and  <j>  ,  supposedly  representing  different  states  of  the 
same  symmetry,  whose  approximate  energy  curves  cross.  In  order  to  get  an 
improved  approximation,  linear  combinations  of  the  original  states  are 
sought  which  diagonalize  the  2x2  electronic  Hamiltonian,  H^.  When  this 
is  done,  the  resulting  new  states  and  are  known  as  "true  adiabatic" 
states.  The  corresponding  true  adiabatic  eigenvalues  E^(R)  and  Ejj(R)  are 
found  to  be  separated  by  the  amount  -  E^  =  |  (H^  “  H ^ ^  +  4Hj^  J 
a  positive  definite  quantity.  In  other  words  they  cannot  cross,  and  must 
behave  in  a  way  illustrated  in  Fig.  1(a).  This  in  brief  is  the  NCR. 

In  addition  to  NCR,  there  is  a  second  and  closely  related  logical 

3 

consequence  of  defining  true  adiabatic  states  which  was  deduced  by  Massey. 
Massey's  theorem  in  essence  points  out  that,  since  by  choice  the  true  adia¬ 
batic  states  diagonalize  H^,  it  follows  that  collision-induced  transitions 
between  electronic  states  cannot  be  mediated  by  a  matrix  element  of  the 
diagonal  operator  Consequently  these  transitions  must  be  caused  by 

T',  the  coupling  kinetic  energy  term  neglected  in  the  Born-Oppenheimer 
approximation,  and  so  this  usn_.~ul  approximation  must  be  abandoned  when 
transitions  occur.  Logical  and  rigorous  as  it  is,  Massey's  theorem,  like 
NCR  in  its  early  days,  has  had  difficulty  of  acceptance  in  t.iat  it  seems 

to  conflict  with  much  of  what  is  known  of  reactions,  especially  the  theory 

4  5 

of  dissociative  recombination  and  attachment.  How  misleading  (though 
formally  correct)  these  two  theorems  are  when  applied  to  collisions  is 


3  H.S.W.  Massey,  Negative  Ions  (Cambridge  Univ.  Press,  N.Y.,  1950). 
A  D.  R.  Bates,  Phys.  Rev.  ]_8,  492(L)  (1950)  . 

5  T.  F.  O'Malley,  Phys.  Rev.  150,  14  (1966). 
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brought  home  by  the  fact  that  even  Massey  was  led  in  1937  to  conclude  on 
their  strength  that  dissociative  recombination  must  necessarily  be  a  very 
slow  reaction.  It  was  not  until  many  years  later^  that  this  was  to  some 
extent  corrected. 

We  shall  next  point  out  very  briefly,  for  those  not  already  acquainted 
with  them,  some  of  the  shortcomings  of  true  adiabatic  states  and  the  con¬ 
sequent  NCR  when  applied  to  collision  problems.  Then  the  alternative  "quasi- 
adiabatic"  states  11  be  introduced,  which  have  the  property  that  they  may 
cross  freely  and  do  not  require  non-adiabatic  (T1)  coupling.  These  states 
should  probably  be  used  in  connection  with  all  problems  involving  collis- 

g 

ions  or  dissociations. 

Many  of  the  problems  and  distressing  features  of  the  NCR  are  illus¬ 
trated  with  the  help  of  Fig. 1(a)  as  contrasted  with  1(b).  The  curves  in  Fig. 
1(b)  mlg*  ;  be  found,  for  example,  when  one  first  does  a  calculation  with  very 
simple  Ueitler-London  or  molecular  orbitrl  functions.  A  slightly  more 
elaborate  calculation,  as  described  above,  would  give  curves  with  the  prop¬ 
erties  of  true  adiabatic  states  as  shown  in  Fig.  1(a)  with  their  avoided 
crossings  and  "obligatory  humps." 

Aside  from  the  obvious  aesthetic  objections  to  the  kind  of  curves  in 
Fig. 1(a),  there  are  a  number  of  more  practical  drawbacks  connected  with  these 
true  adiabatic  curves.  First  there  is  the  fact  if  the  colliding  atoms  h  ve 
almost  any  finite  velocity  the  probability  of  their  actually  following  the 
true  adiabatic  curves  as  opposed  to  the  crossing  curves  is  known  to  become 
vanishingly  small,  as  is  indicated  in  the  next  paragraph.  Next,  the  Born- 
Oppenheimer  approximation,  which  is  generally  so  very  useful  and  accurate 


°  H.S.W.  Massey,  Proc.  Roy.  Soc..  A163.  542  (1937). 

^  D.  R.  Bates  and  H.S.W.  Massey,  Proc.  Roy.  Soc.  A187 .  261  (1947). 

®  Dissociation  is  here  taken  to  mean  the  dissociation  of  repulsive 

states  mainly,  or  predissoci.af.ion  by  curve  crossing.  Predissociation 
by  barrier  tunnelling  is  another  matter. 
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for  molecular  states,  breaks  down  completely  in  the  neighborhood  of  an 
avoided  crossing  when  one  uses  true  adiabatic  states.  This  follows  from 
Massey's  theorem.  Finally,  tlu  -e  are  two  large  classes  of  collisions  where 
theorists  have  long  since  been  forced  to  abandon  the  NCR,  at  least  for  prac¬ 
tical  purposes,  and  to  think  in  terms  of  tie  crossing  curves  in  Fig. 1(b). 

The  first  class  of  collisions  where  the  NCR  has  hid  to  be  suspended 

is  the  slow  reactive  or  inelastic  collisions  between  heavy  particles  which 

9 

are  well  described  in  the  low  energy  region  by  the  Landau-Zener  theory, 
together  with  subsequent  improved  versions.  In  deriving  this  theory,  one 
starts  with  two  "approximate"  electronic  states,  which  are  never  defined 
precisely,  but  whose  physical  interpretation  is  assumed  to  be  clear.  These 
states  have  potential  curves  which  "pseudo-cross."  The  transition  proba¬ 
bility  of  passing  from  one  of  the  pseudocrossing  curves  to  the  other  is 

■"26 

derived,  after  some  approximations,  and  found  to  be  1  -  e  for  each  pas- 
2 

sage,  where  6  =  n  V  /hv  A F.  (V  is  the  energy  of  interaction  between  the 

crossing  states,  so  that  2V  would  be  the  smallest  separation  of  the  true 

adiabatic  states;  v  is  the  internuclear  velocity  and  AF  is  the  difference 

of  slopes.)  One  consequence  of  the  Landau-Zener  formula,  as  mentioned 

“26 

above,  is  that  for  typical  parameters  the  probability  1  -  e  of  follow¬ 
ing  the  true  adiabatic  curve  becomes  negligible  for  almost  any  finite 
velocity.  Now  although  this  theory  has  been  found  to  describe  accurately 
all  sorts  of  low  energy  reactions,  it  is  still  necessary  to  apologize  for 
the  "pseudocrossing"  (which  is  essential  to  the  derivation)  and  to  point 
out  that  the  true  adiabatic  curves  of  course  cannot  really  cross.  Perhaps 
of  more  practical  concern  is  the  fact  that  the  pseudocrossing  electronic 
states  simply  are  not  defined,  except  for  being  "sufficiently  inaccurate" 
to  allow  them  to  cross.  A.s  better  computer  calculations  become  feasible, 
this  shortcoming  should  become  progressively  more  and  more  of  a  problem. 

The  second  class  of  collisions,  which  cannot  even  be  described  con¬ 
ceptually  in  terns  of  the  true  adiabatic  states  and  the  NCR,  consists  of 

^  L.  Landau,  Physik.  Z.  Sovjetunion  2_,  46  (1932).  C.  Zener,  Proc.  Roy. 

Soc.  A137.  696  (1932). 
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the  processes  of  dissociative  recombination  (e  t  AB  AB  -*  A  +  B) ,  and 

5  »  r 

dissociative  attachment  (e  +  AB  .  AB^  -*  A  +  B  ).  (Fortunately,  these  pro¬ 
cesses,  when  prrperiy  analyzed,  also  point  the  way  to  the  needed  definition 
cf  "quasi-adiabatic  states.")  The  recombination  process  is  understood  by 
reference  to  Fig. 1(b).  An  electron  collides  with  the  ion  AB+,  and  causes 
an  electronic  transition  to  the  state  labeled  (AB)^,  which  state  finally 
dissociates  to  A  +  B.  It  is  seen  that  the  curve  (AB)^,  in  coming  in  from 
=  first  had  to  violate  NCR  countless  times  and  cross  an  entire  Rydberg 
series  of  the  same  symmetry  before  reaching  the  curve  AB+.  More  seriously, 
after  crossing  AB+  and  entering  the  hatched  region  (the  electronic  continu¬ 
um  of  AB) ,  the  curve  must  then  cu'-  through  a  continuum  of  states,  in  flag¬ 
rant  violation  of  NCR.  This  latter  situation  is  exactly  the  same  for  the 
closely  related  dissociative  attachment  process. Figure  1(a)  is  the  way 
these  curves  would  have  had  to  be  drawn  had  the  NCR  not  been  ignored. 

Ironically  the  crossing  of  the  continuum,  which  presents  the  strongest 
violation  of  NCR,  also  points  the  way  to  the  solution  of  the  entire  non¬ 
crossing  dilemma.  For  a  discrete  state  like  AB^  which  is  degenerate  with 
(crosses)  a  continuum,  is  a  familiar  situation  in  theoretical  physics  and 
has  been  dealt  with  exhaustively.^  Such  a  discrete  state  is  the  reson¬ 
ance"  or  "quasi-stationary  state."  It  is  by  use  of  such  states  that  the 
problem  of  dissociative  attachment  (and  by  inference  recombination)  has 
been  successfully  treated  in  a  fully  quantitative  way."* 


The  definition  of  the  resonance  or  quasi-stationary  state  for  the 
autoionizing  electronic  states  c.  molecules  is  discussed  in  some  detail  in 

^  For  the  present  purpose,  the  main  features 


Ref.  5  as  well  as  elsewhere 
are  as  follows 

artifically  restricted  in  some  way  to  a  portion  of  Hilbert  space  by  a  pro- 


The  electronic  wave  function  -y  of  the  resonance  state  is 

r 


10 


11 


H.  Feshbach,  inn.  Phys.  (New  York)  5,  357  (1938);  19,  287  (1962);  and 
references  to  previous  work  therein  contained. 

.  N.  B;  rdslcy,  A.  Herzcnberg  and  F.  Mandl.  Ato— lr  Collision  Prnresse-~ 
(M.R.C.  McDowell,  Ed.)  North  Holland,  Amsterdam  (1964),  p.  415.  .I.C.Y, 
Chen,  J.  Chem.  Phys.  40,  3507  (1964). 
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jection  operator  Q^o) ,  while  its  energy  eigenvalue 

<q(0)  VHel!  Q(0)  V  =  C(R)  (1) 

is  made  stationary.  (As  in  any  adiabat  c  definition,  R  appears  parametric¬ 
ally.) 


Having  defined  ^ ,  a  new  projection  operator,  Q  =  4f  ><:  4r  is  con¬ 
structed  which  together  with  P  -  1  -  Q  decomposes  the  electronic  Hamiltonian 

H  1  as  follows 
el 


H 


el 


PH  ,P 
el 


PHeiQ 


QHeiP  QHelQ 


(2) 


The  eigenstate  corresponding  to  Pi.^P,  called  4p  the  "potential  scattering" 
state,  is  constrained  by  P  to  be  orthogonal  to  4  ■  The  above  decomposition 
of  H^^ ,  or  more  properly  of  the  full  H,  f  rmed  the  basis  of  the  new  formal 
theory  of  rearrangement  collisions  introduced  in  Ref.  5. 


The  two  states  4r  and  4p  have  been  defined  to  lie  in  orthogonal  por¬ 
tions  of  Hilbert  space  and  they  leave  Hg^  undiagonalized  (note  that  it  is 
only  the  full  H  which  must  be  diagonal).  One  consequence  of  a  non-diagonal 
H^  is  that  transitions  from  4^  to  4^_  may  be  caused  by  the  off-diagonal 
electronic  matrix  element  PH^Q,  ir  contrast  to  the  Massey  theorem.  Thus 
the  coupling  term  T'  is  no  more  important  here  than  in  ordinary  molecular 
states  and  the  usefulness  of  the  Born-Oppenhcimer  approximation  is  preserved. 
More  important,  end  in  fact  paramount  in  the  present  treatment,  the  mathe- 


combined  as  were  the  states  of  von  Neumann  and  Wigner,  into  "true  adiabatic" 
states  of  H^,  resulting  in  the  NCR.  Rather  they  are  allowed  to  serve  as 
a  convenient  starting  point  for  an  exact  solution  cl  the  collision  problem. 
(This  theory  will  be  published  elsewhere.)  By  thus  preserving  their  ident¬ 
ity,  the  energies  bp  and  of  these  states  are  not  constrained  to  avoid 
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crossing,  and  the  levels  can  and  do  err  .  The  lack  of  a  prohibition 
against  crossing  may  perhaps  be  made  clear  by  noting  that  as  eigenvalues 
of  different  Hamiltonians,  PHP  and  QHQ,  the  two  levels  are  entirely  inde¬ 
pendent  of  one  another. 

What  has  been  said  so  far  of  the  quasi-stationary  states  and  d>^ 
has  referred  explicitly  to  the  hatched  region  of  Fig. 1(b),  i.e.,  to  the 
region  where  the  state  AB^,  with  wave  function  4>r»  is  autoionizing.  It 
is  the  burden  of  this  paper  to  point  out  that  all  the  difficulties  result¬ 
ing  from  true  adiabatic  states,  i.e.,  the  NCR  and  the  Massey  theorem,  are 
avoided  if  only  one  uses  exactly  the  same  definition  of  states  (with  exactly 
the  same  projection  operators  and  Q)  for  all  R  even  after  the  state 

AB  has  fallen  below  the  continuum  and  can  no  longer  autoionize.  It  is  not 
r 

necessary  to  make  any  changes  at  all  in  the  discrete  region.  Again  the 

(o' 

quasi-stationary  state  is  defined  with  the  same  projection  operator  Q  '  . 
The  resulting  new  projection  operator  P  is  then  used  to  define  the  ortho¬ 
gonal  "potential  scattering"  function  Again  the  corresponding  potential 

curves  and  are  not  constrained  by  the  NCR  or  by  Massey's  theorem. 

The  only  minor  difference  is  that  the  name  "potential  scattering"  function 
is  no  longer  very  appropriate,  since  ^  is  now  a  discrete  quasistationary 
state,  on  the  same  footing  with  $  . 


The  name  which  is  suggested  for  the  crossing  states  and  $  ,  defined 
as  above  and  extended  to  the  discrete  region,  is  quasi-adiabatic  states, 
short  for  quasi-stationary  adiabatic  states.  This  is  appropriate  since  they 
are  adiabatic,  being  defined  w'.th  R  fixed,  but  unlike  the  "true  adiabatic" 
states  they  are  only  quasi  or  artificial  bound  states. 

To  avoid  possible  misunderstanding,  something  should  be  said  here 

about  what  is  presently  meant  and  not  meant  by  the  projection  operator  . 

The  obvious  definition  of  this  projection  operator  and  presumably  the  most 

i  o 

useful  Is  that  of  Feshbach,  which  defines  a  state  in  which  an  electron  is 
bound  to  an  excited  core.  However,  is  meant  to  he  somewhat  more  flex¬ 

ible.  For  one  thing,  certain  simplifications  <f  F..ubl'ach's  states  may  be 
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1  2 

desirable.  In  general,  any  artificial  restraint  on  the  wave  function 
which  leads  to  a  physically  neaningful  and  useful  state  is  considered  to  be 
included  as  a  possible  since  this  tiny  always  be  thought  of  as  being 

produced  by  some  kind  of  projection  operator.  As  one  example,  an  artificial 
constraint  might  sometimes  be  useful  for  defining  ionic  states.  However, 
at  present,  only  the  Feshbach  type  of  core  excited  state  is  actually  en¬ 
visaged. 

There  are  a  number  of  examples  in  the  literature  of  states  which  have 

been  defined  in  such  a  way  that  they  happen  to  satisfy  the  definition  of 

quasi-adiabatic  states.  Two  of  these  will  be  mentioned.  First,  there 
2  1 

is  the  (o  Is)  E  state  of  H„,  as  calculated  with  the  simplest  single  con- 
u  8  13  ^  2 

figuration  wave  function.  This  configuration  (culs)  satisfies  the  defi¬ 
nition  of  the  lowest  Feshbach-type  resonance  in  the  H0  system,  and  so  is  a 
quasi-adiabatic  state.  The  potential  curve  of  thx  state  has  been  plotted 
as  an  illustration  in  Ref.  12.  As  expected,  it  cuts  right  across  higher 
states  of  the  same  symmetry,  crosses  the  entire  Rydberg  series  and  ulti¬ 
mately  becomes  autoionizing,  just  as  AB^  of  Fig  1(b).  An  interesting  inci¬ 
dental  feature  of  the  curve  is  that  it  becomes  ionic  at  large  R.  Naturallj , 
as  more  ambitious  calculations  have  been  made,  looking  for  true  adiabatic 
states,  this  state  has  been  mixed  with  others,  resulting  In  humps  and 
avoided  crossings  as  in  Fig.  1(a). 


A  second  example  of  quasi-adiabatic  electronic  states  in  the  litera¬ 
ls 

ture  is  seen  in  the  He,,  states  found  by  Mulliken  and  labeled  B  Core  states. 

“  3  +  3  + 

These  included  the  B2s  E ,  and  a  B3p  £u  state.  The  concept  of  B  Core 
state,  as  pointed  out  in  Raf .  4,  directly  satisfies  Feshbach’s  definition 
of  a  resonance  state.  Consequently  we  expect  first  that  these  states  should 


^  See  for  example  L.  Lipsky  and  A.  Russek,  Phys.  Rev.  142,  39  (1966). 

^  W.  Kauzmann,  Quantum  Chemistry  (Academic  Press,  New  York,  1957),  p.  397. 

14  -  .  ...  -  .  ..  ..... 

~  '  ft*  a.  iiuiiiKeu,  ruyt».  t\ev*  . 
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be  autoionizing  for  small  R,  and  in  fact  they  are.  Further,  as  R  increases 

they  should  cut  through  the  Rydberg  series  of  the  same  symmetry,  which  they 

also  do.  However,  both  of  these  curves  have  been  drawn  in  Ref.  14  so  as  to 

3  +  3  + 

avoid  crossing  the  lowest  state  of  the  series  (  Z  and  £  )  and  so  they 

o  g  u 

both  turn  up  sharply  above  R  =  2A,  somewhat  in  the  fashion  of  Fig. 1(a).  It 
is  suggested  here  that  if  the  definition  of  quasi-adiabatic  state,  implied 
by  the  B  Core  definition,  is  applied  consistently,  then  these  states  need 
not  avoid  crossing  even  the  'owest  states  in  the  series,  but  it  is  possible 
that  they  would  simply  decrease  monatonically  with  R  as  the  AB  curve  of 
Fig.  1(b),  and  dissociate  into  the  lowest  S  and  P  state  of  lie  respectively. 


This  last  suggestion  would  tend  to  remove  the  dis  irepancy  which 

presently  exists  between  these  curves  of  Mulliken  and  subsequent  results 

of  Robertson. ^  The  latter  observed  a  state  of  the  system  which  is  auto- 

3 

ionizing  at  small  R  and  which  dissociates  into  the  lowest  P  state  of  He 
in  apparent  disagreement  with  Mulliken* s  curves.  But  if  the  definition  of 
quasi-adiabatic  B  core  states  is  adopted  as  is  here  suggested,  then  for 
any  given  states  of  the  dissociation  products  (and  in  particular  for  those 
found  by  Robertson)  there  ought  to  be  a  corresponding  molecular  curve  satis¬ 
fying  the  definition  of  a  B  core  quasi-adiabatic  state.  This  curve  would 
he  expected  to  be  purely  repulsive,  following  the  shape  of  the  B  core,  and 
at  smaller  R  to  cross  the  Rydberg  series  and  become  autoionizing.  Thus, 
the  existence  of  a  suitable  potential  curve,  corresponding  to  Robinson's 
data,  should  present  no  problem  once  the  crtificiai  restraints  imposed  by 

NCR  arc-  removed  by  using  quasi-adiabatic  states.  It  is  not  clear  at  present 

3 

whether  this  particular  state  might  be  the  E  curve  of  Mulliken,  redrawn 

3^  3  u 

to  join  B  Z  ,  or  perhaps  a  li  B  core  state,  similarly  redrawn. 


In  summary,  a  definition  has  been  given  of  quasi-adiabatic  states 
which  It  is  proposed  should  replace  the  conventional  true  adiabatic  states 

Q 

for  all  processes  involving  collisions  or  dissociations.  In  particular, 
it  was  shown  that  such  states  have  the  advantage  of  not  being  burdened  by 
a  non-crossing  rule,  so  that  they  may  reproduce  the  curves  that  are  actually 
followed  in  a  physical  collision  or  dissociation.  Fm-tln-r  no  not 


w.  W.  Robertson,  .1.  Chen,.  Phys.  4_2,  2064  (1965);  C.  li.  Collins  and 
W.  W.  Robertson,  J.  Chem.  Phys,  4_3,  4188  (1965). 
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satisfy  Massey's  theorem,  so  that  the  Born-Oppenheimer  approximation  may 
be  retained.  These  states  may  be  of  most  immediate  usefulness  in  connection 
with  the  semi-empirical  process  of  constructing  potential  curves  to  fit  or 
explain  given  experimental  data.  An  improved  theory  of  heavy  particle  col¬ 
lisions  based  on  these  states  will  be  published  elsewhere. 

Finally,  these  states  are  not  meant  to  replace  the  conventional  true 
adiabatic  states  in  those  areas  where  the  latter  are  actually  useful,  namely 
where  such  states  are  either  very  widely  separated  or  parallel  and  where  in 
addition  the  discrete  vibrational  level  of  interest  is  sufficiently  low  so 
that  the  Landau-Zener  exponent  6  would  be  very  large. 

Thanks  are  due  B.  A.  Lippmann  and  P.  J.  Redmond  for  stimulating  dis¬ 
cussions. 
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THEORY  OF  RESONANT  DISSOCIATIVE  RECOMBINATION 

by 

Peter  J.  Redmond 


The  theory  is  developed  with  particular  emphasis  on 
the  interaction  of  the  resonant  state  and  the  infi¬ 
nite  set  of  Rydberg  states  crossed  by  the  resonant 
state.  It  is  shown  that  the  interaction  is  large 
only  for  a  finite  subset  of  the  Rydberg  states. 


In  this  paper  we  consider  the  process  e  +  AB  -*■  A  +  B  where  the 
atoms  in  the  final  state  will  generally  be  excited.  We  shall  show  that 
the  total  cross  section  corresponding  to  a  sum  over  all  possible  exit 
channels  can  be  calculated  in  a  fairly  straightforward  way.  In  principle 
the  theory  also  permits  the  calculation  of  branching  ratios  although  such 
a  calculation  would  be  much  more  involved. 

We  assume  that  the  reaction  proceeds  by  the  formation  of  a  resonant 

*  +  .  * 

complex  AB  so  that  e  +  AB  Z  (AB)  -*■  A  +  B.  It  is  well  known  that  there 
are  an  infinity  of  such  resonant  states  lying  below  the  first  electronic¬ 
ally  excited  state  of  the  AB+  system.  It  is  also  known  that  the  widths  of 
the  resonant  states  are  small  compared  to  the  spacing  of  the  resonant  states 
so  that  they  are  excited  independently.  However,  when  the  atoms  A  and  B 
separate  they  move  along  a  potential  curve  which  must  cross  (or  "pseudo 
cross")  an  infinite  Rydberg  series  of  states.  These  curve  crossings  may 
lead  to  a  coupling  between  the  resonant  states  and  will  lead  to  a  distri¬ 
bution  of  the  particles  over  a  variety  of  exit  channels.  It  is  this  latter 
phenomena  which  makes  the  theory  of  dissociative  recombination  more  compli¬ 
cated  than  the  theory  of  dissociative  attachment.  We  shall  show  that  only 


Submitted  for  presentation  at  the  Fifth  International  Conference  on 
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a  finite  (and  usually  small)  number  of  the  resonant  states  can  bo  coupled 
ln  this  fashion  and  that  the  coupling  can  be  described  by  an  energy  depend- 
ent  matrix  for  a  reflection  amplitude. 

The  wave  function  for  the  system  may  be  written  in  the  form 

^(R,P*k)  =  Zy  AO,R)  <J(*»R»p) 


+  B  (R)  x(+)  (R*P.k) 

o\  o  % 


+  fd k'  l.  B(j,k‘,R)  x(  ^(J.k'.R.o) 

%  j  A.  -V  ^  ^ 


(1) 


where  R  is  the  nuclear  separation  and  e  represents  the  electronic  coordi¬ 
nates.  The  wave  functions  H')  are  the  electronic  wave  functions  for  the 
resonant  states.1  The  X  are  the  electronic  scattering  states  and 
it  is  assumed  that  they  diagonalize  the  electronic  part  of  the  Hamiltonian 
in  the  subspace  of  states  orthogonal  to  the  ♦(A).  The  X1  satisfy  incom¬ 
ing  boundary  conditions  and  the  x“  outgoing  conditions  and  the  vector  k 
represents  the  momentum  of  the  incoming  or  outgoing  electron.  The  term 
B  x(+)  is  the  incoming  wave  and  Bq  is  the  nuclear  wave  function  of  the 
initial  AB+  state.  The  states  ♦(A)  and  X  are  slowly  varying  functions  of 
the  nuclear  separation.  This  is  possible  only  if  the  energy  curves  E, (R) 
corresponding  to  the  $(*)  are  allowed  to  cross. 

The  only  non-diagonal  matrix  elements  of  the  electronic  Hamiltonian 
are  then 


<  *(X,R),  Hel,  X(+)(j,k',R)  >  =  Y  (A, j .k' ,R) 


(2) 


and 


<  $(\,R),  H  L ,  f(u,R)  ’  W(a,u,R) 


(3) 


1 


We  use  the  term  resonant  state  to  include  all 
*- r» q  rtpufr.il  qvfifeir.  which  interact  with  the 

\J  t  v  *  ■  •  y 


those  discrete  states 
state  initially  excited. 
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The  Y  and  W  are  non-vanishing  in  separated  regions  of  configuration 
space.  The  matrix  element  W  is  appreciable  only  in  the  immediate  vicinity 
of  a  curve  crossing  whereas  the  y's  are  non-vanishing  only  when  the  reson¬ 
ant  state  curve  lies  above  the  ground  state  potential  curve  of  AB  . 

Since  the  $  and  x  a re  slowly  varying  functions  of  R,  the  many  state 
Bom-Oppenheimer  (adiabatic)  approximation  is  appropriate  and  (E-H)  <J>  =  0 
is  equivalent  to  the  set  of  coupled  equations. 

(E  -  E,  (R)  +  ~  V2)  A(A,R)  =  Y, (A,0,k,R)  B  (R) 

+  /  dk'  Z.  B(j,k',R)  Y  (A.j.k'.R) 

j  'v  'v  “ 

+  Z  W(A,u,R)  A(u »R)  (4) 

Va  f\j 

and 

2 

(e  -  E^R)  -  ^  v2)  B(j,k’,R)  =  Ex  A(A,R)  Y*(A,j,k',R)  , 

(5) 

where  m  and  M  are  the  electronic  and  reduced  atomic  masses  respectively. 

The  solution  of  Eq.  (5)  which  satisfies  the  proper  boundary  conditions 
may  be  written  asymptotically 


B(j,k',R) 


A(A,R)  Y*(A,j,k') 

^  'V 


(6) 


If  the  Green's  function  in  Eq.  (6)  is  expanded  in  terms  of  the  eigenstates 
of  the  nuclear  Hamiltonian,  the  discrete  states  describe  elastic  scatter¬ 
ing  and  vibrational  ana  rotational  excitation.  The  continuum  states  con¬ 
tribute  to  dissociation  although  the  main  contribution  to  dissociation 
comes  irom  the  ACA.R)  terms.  For  example,  the  elastic  scattering  comes 

'v 

from  the  term 
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b  w  „ 

O  j  .0 


B  (j ,k* ,R)  «  — - -  „ 

elastic  %  ^  ^2  ^,2 

2m  2m 


£.  I  dR'  B  (R ' )  A.(R')  y. (A.O.k'.K'). 

■  O  X  A  -V  +  1, 


S-  it: 


(7) 


(+1 

This  term  interferes  with  the  outgoing  wave  part  of  B  (R)  x  ' (0,k,R)  to 

o  x 

give  a  typical  resonance  pattern  for  the  elastic  scattering. 

In  <  to  determine  the  dissociative  cross  section  we  take  the 

divergence  or  the  heavy  particle  flux.  We  find 


I.(A*(>.R)  7  A(X,R)  -  A(X,R)  7A*(XfR)) 

A  'v  x  -c  -c 


+  Z  I  dk’ (B*(j,k’ ,R)  7  B  (  j  ,  k’ , R)  -  B(j,k',R)  V  B*(j,k',R))  I 

J  V  XT.  \  %  XX  'V'Vl 


-  1  (A  (A,R)  Y . ( X , 0 , k , R )  B  (R)  -  A(A,R)  y  (A.O.k.R)  B  (R))  .  (8) 

A  x  +  O  %  -v  +  O  v 


Since  there  is  no  incoming  heavy  particle  flux  the  total  cross  section 

for  dissociation  can  be  obtained  by  integrating  the  right-hand  side  of 

Eq.  (8)  over  all  space.  Because  of  the  presence  of  the  factor  B  (R)  on 

o 

the  right-hand  side  of  Eq.  (8)  (and  in  the  integral  occurring  in  Eq.  (7)), 
it  is  only  necessary  to  determine  the  functions  A(X,R)  in  the  Franck- 
Condon  region. 

In  determining  A(X,R)  it  is  necessary  to  consider  two  different 

x 

possibilities.  If  the  electron  energy  is  sufficiently  high  so  that  all 
the  attractive  potential  curves  which  are  crossed  by  the  resonant  curve 
involved  correspond  to  exit  channels  which  are  open,  then  any  outgoing 
flux  generated  on  the  Franck-Condon  region  will  remain  outgoing.  The 


2  T.  F.  O’Malley,  Phys.  Rev.  150.  14,  1966. 

^  J.  N.  Bardsley,  A.  Herzenberg,  and  F.  Mandl,  Atomic  Collision 

Processes .  edited  by  M.R.C.  McDowell  (North  Holland  Publishing  Co., 
Amsterdam  i964),  p.  145. 
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A(A,R)  can  then  be  determined  in  precisely  the  name  way  that  O'Malley  ^ 

3 

and  Mandl  and  Herzberg  treated  the  dissociative  attachment  problem. 

The  only  difference  is  that  the  final  formula  in  the  recombination  case 
corresponds  to  several  exit  channels  instead  of  only  one  as  in  the  attach¬ 
ment  problem. 

Ii.  some  of  the  exit  channels  with  attractive  potential  curves  are 
closed,  some  of  the  outgoing  flux  generated  in  the  Franck-Condon  region 
will  correspond  to  particles  which  are  temporarily  bound  and  return  to 
the  Franck-Condon  region.  The  finite  number  of  resonant  states  coupled 
to  such  exit  channels  will  then  be  coupled  to  each  other.  These  effects 
can  be  described  by  a  finite  reflection  coefficient  matrix  which  gives 
the  amplitude  of  the  reflected  wave  in  any  resonant  state  produced  by  an 
outgoing  wave  in  another  coupled  resonant  state.  T..^a  matrix  can  be 
determined  by  using  the  Landau-Zener  theory  (or  some  variant  of  it)  and 
considering  the  possible  histories  of  outgoing  particles. 

The  system  e  +  H ^  •*  H  +  H  provides  a  comparatively  simple  illustra¬ 
tion  of  the  various  possibilities  discussed  above  and  the  theory  will  be 
applied  to  the  system  at  the  meeting. 
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FORMULATION  OF  THE  TWO-STATE  PROBLEM 
IN  TERMS  OF  ADIABATIC  POTENTIAL  CURVES 

by 

Peter  J.  Redmond 

ABSTRACT 


When  two  molecular  adiabatic  potential  curves  nearly 
cross,  the  electronic  wave  function  is  a  rapidly  vary¬ 
ing  function  of  the  intemuclear  separation.  The  one- 
state  Born-Oppenheimer  approxiriation  is  then  not  valid. 

By  introducing  appropriate  linear  combinations  of  the 
two  adiabatic  states  involved  a  two-state  Bom-Oppen- 
heimer  approximation  can  be  made .  The  resulting  pair 
of  coupled  equations  for  the  nuclear  wave  function  are 
identical  with  those  obtained  in  the  conventional  two- 
state  theories.  However,  all  the  parameters  can  be 
determined  by  an  inspection  of  the  adiabatic  curves 
and  the  validity  of  the  approximation  scheme  can  be 
justified. 

I.  INTRODUCTION 

In  the  usual  approach  to  molecular  problems,  or  to  collision  problems, 
one  makes  the  adiabatic  (Born-Oppenheimer)^  approximation.  The  Hamiltonian 
is  written  as  the  sum  of  a  nuclear  kinetic  energy  term  and  an  electronic 
Hamiltonian 


H 


1_ 

2M 


(R) 


(1) 


and  the  eigenstates  of  H^  are  calculated 

H  4>  (R.O  =  E  (R)  <fc  vR.O  .  (2) 

ei  a  a  a 


*  Submitted  to  ARPA  for  publication  clearance. 

*  M.  Born  and  J.  R.  Oppenheimer,  Ann  der  Phys  84.  457  (1927). 
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The  total  wave  function  may  be  written  as  a  linear  combination  of  these 
electronic  states  with  coefficients  which  are  functions  of  the  internuclear 

separation  R  , 

% 

*  -  E  F  (R)  *  (K.O  •  (3) 

&  01  01 

The  equation  for  F  (R)  is  then 

a  % 

(-k  ’2  +  Ea<R)  *  E)  Fa  ' 

Eb(h  <  ♦«<*>•  ”b(R)  ’  •  ',F6(R)  +  2M  '  ♦»<*>•  ,2  VR?  ’  F6(R))- 

(4) 

The  Born-Oppenheimer  approximation  consists  of  setting  the  right-hand  side 

of  Eq.  (4)  equal  to  zero  (although  the  diagonal  terms  may  fruitfully  be  re- 
2 

tained) .  The  equations  for  the  various  F^  then  decouple.  We  shall  •  .r 
to  this  situation  as  the  one-state  Born-Oppenheimber  approximation.  iuse 
of  the  largeness  of  the  atomic  to  electronic  mass  ratio  this  is  a  good  ap¬ 
proximation  unless  the  electronic  wave  functions  <|>a  exhibit  an  anomalously 
rapid  variation  as  the  internuclear  separation  is  changed.  We  shall  show 
that  this  is  the  case  when  two  of  the  adiabatic  curves  with  the  same  sym¬ 
metry  come  very  close  to  one  another  ("near  crossing"). 

When  there  is  a  near  crossing  it  is  generally  argued  that  the  adia¬ 
batic  states  no  longer  provide  a  suitable  basis  for  an  expansion  of  the 

wavefunction.  The  problem  is  then  reformulated  either  in  a  time  dependent 

3  4 

way  using  atomic  wave  functions  or  in  a  time  independent  way  using  simple 

2  L.  D.  Landau  and  E.  M.  Lifshitz,  Quantum  Mechanics  (Pergamon  Press, 
London-Paris,  1958)  p.  275. 

3  C.  Zener,  Proc.  Roy.  Soc.  (London),  A137.  696  (1932);  L,  Landau, 

Phys.  Z.  Sowjet  Union,  2_,  46  (1932);  D.  R.  Bates,  H.  C.  Johnston 
and  I.  Stewart,  Proc.  Phys.  Soc.  (London)  84,  517  (1964);  C.  A. 

Coulson  and  K.  Zalewski,  Proc.  Roy.  Soc.  (London)  A268,  437  (1962); 

D.  R.  Bates,  Proc.  Roy.  Soc.  (London)  A257.  22  (1960). 

^  N.  F.  Mott  and  H.S.W.  Massey,  The  Theory  of  Atomic  Collisions.  Third 
Ediction  (Oxford  University  Press,  London,  1965)  Chv.p.  13. 
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molecular  orbitals. In  both  of  these  formulations  the  resulting  potential 
curves  cross  and  a  set  of  coupled  equations  are  considered  which  involve 
the  two  states  that  cross. 

There  are  two  disadvantapas  to  this  approach.  Because  of  the  crude 
wave  functions  used  the  electronic  Hamiltonian  will  have  large  matrix 
elements  connecting  different  states  and  it  in  therefore  difficult  to 
justify  the  use  of  only  two  electronic  states  m  formulating  the  problem. 

In  addition  there  is  the  esthetic  disadvantage  that  the  basis  of  the  theory 
disappears  if  the  wave  functions  are  sufficiently  improved  so  that  the 
curves  no  longer  cross. 

In  this  paper  we  reformulate  the  problem  starting  from  the  adiabatic 
potential  curves.  The  resulting  equations  for  the  nuclear  wave  functions 
ave  identical  to  those  obtained  in  the  conventional  approach.  However, 
all  the  parameters  for  the  theory  can  be  obtained  from  an  inspection  of  the 
adiabatic  curve.  In  addition,  it  is  easy  to  justify  the  use  of  a  two-state 
approximation  when  such  an  approximation  is  applicable. 

II.  THEORY 

The  general  shape  of  the  adiabatic  energy  curves  for  a  near-crossing 

is  illustrated  by  the  solid  curves  labeled  Eq  and  Eg  in  Fig.  1.  The  general 

feature  of  such  potential  curves  is  that  the  second  derivative  of  the  energy 

2  2 

with  respect  to  the  iuternuclear  separation  (d  EQ/dR  )  is  large  in  the 
vicinity  of  the  near  crossing  and  that  the  curves  asymptotically  approach 
regions  where  the  second  derivative  is  much  smaller.  We  shall  show  that  in 
this  transition  region  the  adiabatic  electronic  wave  functions 
and  CR* are  rapidly  varying  functions  of  the  parameter  R  with  the  two 


5  W.  Lichten  in  IV th  International  Conference  on  the  Physics  of.._Elg.C^. 
tronic  and  Atomic  Collisions  (Science  Bookcrafters  Inc.,  Hastings  on 
Hudson,  1965)  p.  63;  V.  Fano  and  W.  Lichten,  Phys.  Rev.  Letters  14, 
627  (1965).  The  molecular  orbital  approach  has  the  great  advantage 
that  it  provides  a  large  amount  of  qualitative  information  about  the 
general  nature  of  the  potential  curves  at  the  cost  of  relatively 
little  effort. 
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wave  functions  essentially  changing  roles  as  the  transition  region  is  cros¬ 
sed. 


The  difference  of  the  slopes  of  the  two  energy  curves  changes  sign 
as  the  transition  region  is  crossed.  The  slopes  must  therefore  be  equal 
at  some  intermediate  point  R*  given  by 

3E  3E. 

if  *  if  '  x  “  R  -  **  •  <5) 

In  the  transition  regie  i  the  electronic  wave  function  $  (R,£  , )  can 

a  ex 

be  approximately  represented  as  a  linear  combination  of  the  wave  functions 
<j.o(R*,£el)  and  <t>g(R*»5el)  s°  Chat6 


$  (R)  «  A<|>  (R*)  +  B<J>_(R*)  . 

o  a  p 

The  secular  equations  have  the  form 


(E  (R*)  +  Xp  -  E  (R))  A  +  uo  B  -  0 
a  a 


and 

UP 

A  +  (Eg(R*) 

+  Xp 

-  E 

c 

t(R)>  B  =  0 

where 

P 

-  R  -  R* 

and 

u 

-  <  4>a(R*) 

|r  • 

3Hi 

3R 

4>g(R*)  >  • 

(6) 


(7) 


(8) 


It  is  a  trivial  matter  to  solve  these  equations  and  we  find 
Eq(R)  »  I  +  Xp  +  ( +  U2P2  ) 

Eg(R)  =  E  +  Xp  -  +  v2o2  ) 


&  In  the  following  we  shall  generally  omit  the  electronic  coordinates  £. 
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<|>  (R)  -  sin  0  <j>  (R*)  +  cos  0  <t>„(R*) 

Ta  a  p 

4>g(R)  =  -  cos  0  $a(RA)  +  sin  0  $g(R*) 


with 

2E  -  E  (R*)  +  Efl(R*)  ,  AE  -  E  (R*)  -  E  (R*) 
a  P  a  p 


and 


The  equations  for  $a(R)  a°d  are  readily  inverted  to  give 

4>a(R*)  =  Sin  0  <J>o(R)  -  cos  e  4>g(R) 

4ig(R*)  -  cos  0  4>aCR)  +  Sin  0  <|>g(R)  • 

The  linear  combinations  of  *a(R)  and  4>g(R)  occurring  on  the  right-hand  side 
of  Eq.  (12)  are  therefore  slowly  varying  functions  of  R.  The  rapid  changes 
in  <f>a(R)  and  4>g(R)  are  compensated  by  the  rapid  variation  of  0(d0/dR  = 
p/A £  at  R  ®  R*)  The  particular  linear  combinations  g^ver.  by  Eq.  (12)  are 
not  the  most  convenient  set  since  they  do  not  asymptotically  approach  the 
adiabatic  states.  The  sum  and  the  difference  of  the  right-hand  sides  of 
Eq.  (12)  provide  more  suitable  combinations.  We  therefore  define 

4>1(R)  =  cos  x  $a(R)  +  sin  X  $g(R> 

4>2(R)  =  -  sin  x  <t>a(R)  +  cos  x  <t>B(R)  ^) 

where  x  ■  6  -  45°.  To  the  left  of  the  transition  region  x  4  0  so  that 

$  -*•  A  and  At  R  -  R*,  X  -  45°.  To  the  right  of  the  transition 

a  1  d  2 

region  X  90°  so  that  <t>B  -  ♦1  and  <J>a  -*2-  Since  ^  and  4>2  do  not  change 
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appreciably  the  functions  4>a  and  have  suffered  large  changes  and  have 

essentially  interchanged.  The  functions  ♦.(R)  and  (R)  provide  a  suitable 

1  1  7 

basis  to  be  used  when  making  the  Born-Oppenheimer  approximation.  The  wave 
functions  ^  and  4>2  can  be  used  to  define  new  potential  curves  indicated 
by  the  dotted  line  in  Fig.  1.  The  diagonal  elements  of  determine  Ej 
and  E2  and  the  off  diagonal  element  determines  E ^  the  interaction  between 
the  two  curves  E.  We  find 

EX(R)  *  cos2x  Ea(R)  +  sin2X  E@(R) 

E2(R)  -  sin2X  Ea(R)  +  cos2X  Eg(R) 

E12(R)  -  +  cos  x  sin  x  (Eg(R)  -  Ea(R))  .  (14) 

We  can  now  make  the  Born-Oppenheimer  approximation  with  confidence 
and  we  get  the  coupled  equations  for  the  nuclear  wave  functions 

(-  "„I  la  r2  k  +  +  Ei(W)  Fi(R)  +  Ei2(R)  f*<r)  "  EFi<R)  ’ 

IR  /MR 

('  Zj  ^  R2  ^  +  +  E2(R>)  F2<R>  +  E12(R)  E1<R)  -  EF2<R>  • 

/MR  /MR 

(15) 

These  are  the  usual  equations  for  the  two-state  problem,  except  that 
in  our  case  the  potentials  are  determined  from  the  adiabatic  curves.  The 
usual  treatments  of  the.  two-state  problem  are  therefore  applicable. 

All  of  the  parameters  needed  to  determine  the  nuclear  wave  functions 
can  be  obtained  directly  from  the  adiabatic  potential  curves  except  for 

^  The  electronic  wave  functions  vary  rapidly  with  internuclear  separa¬ 
tion  in  another  circumstance  -  when  there  is  an  electronic  resonance. 

The  Born-Oppenheimer  approximation  can  then  not  be  directly  applied 
and  the  appropriate  treatment  is  quite  analogous  to  ours  with  the 
resonant  state  playing  a  role  similar  to  and  See  T.F.  O’Malley, 

Phys.  Rev.  150,  14  (1966).  1  L 
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the  parameter  which  is  defined  hy  the  matrix  element  in  Eq.  (8).  However, 
if  the  near  curve  crossing  is  sufficiently  close  so  as  to  invalidate  the 
use  of  the  one  state  Born-Oppenheimer  approximation,  y  can  be  obtained  to 
sufficient  accuracy  by  fitting  the  energy  curves  in  the  transition  region 
with  the  hyperbola 

(E  (R)  -  Ea (R) ) 2  -  (E  (R*)  -  E.(R*))2  =  Ay2  (R  -  R*)2  .  (16) 

a  d  a  B 
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TREATMENT  OF  ANGULAR  MOMENTUM  IN  MOMENTUM  SPACE* 

by 

Peter  J.  Redmond 

ABSTRACT 

Equations  are  derived,  in  a  momentum  space  repre¬ 
sentation,  for  the  wave  function  of  a  many  body 
system  with  a  definite  angular  momentum.  Invariance 
of  the  Hamiltonian  under  time  reversal,  parity,  and 
exchange  has  simple  consequences  for  a  three-body 
system. 


INTRODUCTION 

The  general  form  of  the  wave  function  of  a  system  having  a  definite 
angular  momentum  has  been  given  by  Wigner.  He  showed  that  the  dependence 
of  the  wave  function  on  the  orientation  could  be  Indicated  explicitly  so 
that  the  dynamical  problem  involved  a  set  of  functions  f^  which  depend 
only  on  quantities  which  are  invariant  under  rotations. 

Although  the  general  form  of  the  wave  function  is  known,  considerable 
labor  is  necessrry  before  equations  relating  the  f^'s  can  be  obtained. 

All  the  treatments  familiar  to  the  author  consider  the  problem  in  config¬ 
uration  space.  The  kinetic  energy  is  then  represented  by  differential 
operators  and  the  central  problem  is  to  express  these  deviations  in  terms 
of  derivatives  with  respect  to  angles  and  derivatives  with  respect  to 
rotational  invariants.  It  is  in  the  nature  of  this  process  of  changing 
variables  that  a  successful  completion  of  the  program  for  a  three-body 
system  is  of  little  assistance  if  one  is  confronted  by  a  four-body  problem. 

It  is  often  convenient  to  consider  a  problem  in  momentum  apace  rather 
than  configuration  space  (for  example,  in  a  scattering  problem) .  it  is  the 


*  Submitted  to  ARPA  for  publication  _lerrance. 
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purpose  of  this  paper  to  derive  tire  equations  describing  a  many  body  system 
with  a  definite  angular  momentum  ±n  the  center  of  mass  frame.  We  believe 
that  the  results  are  structurally  much  simpler  than  those  obtained  in  con- 
.ation  space.  In  particular,  tht  equations  have  the  same  structure 
for  a  four-or-more-body  problem  as  they  do  for  a  threv-body  problem. 

The  paper  is  divided  into  four  sections.  Section  1  is  the  introduc¬ 
tion  which  you  have  almost  finished  reading.  In  Section  2  a  thumbnail 
sketch  of  some  aspects  of  the  general  theory  of  angular  momentum  is  pre¬ 
sented  in  order  to  establish  a  notation  and  a  set  of  conventions.  The 
third  section  is  devoted  to  a  derivation  or  our  central  result.  In  the 
last  section  the  effects  of  parity,  exchange,  and  time  reversal  invariance 
are  discussed  for  a  three-body  system. 

SOME  ASPECTS  OF  THE  GENERAL  THEORY  OF  ANGULAR  MOMENTUM 

The  effect  i  of  rotations  may  be  described  from  two  different  points 
of  view.  Acccrd  ng  to  the  first  point  of  view  one  considers  the  change 
in  description  of  a  fixed  geometrical  configuration  produced  by  a  rotation 
of  the  coordinate  axes.  According  to  the  second  point  of  view  rotations 
are  used  to  describe  the  motion  of  a  geometrical  configuration  relative  to 
a  fixed  set  of  coordinate  axes.  Throughout  this  paper  we  adopt  the  second 
point  of  view. 

We  consider  fu-  :tions  ip ( { p } )  of  the  set  of  vectors  p^,  P£»  —  P^ 

(*  {p})  which  we  denote  by  an  inner  product 

ij/({p})  =  <  ip }  | ip  >  (1) 

An  ii.f initesmal  rotation  about  the  r.  axis  through  an  angle  6<j>  causes 
a  change  in  the  vector  p^  by  an  amount  6p^  given  by 

6p  -  6<p  n  x  p  (2) 

a  ra 

This  causes  an  inf  initesmal  change  in  the  function  iKtp-)  which  muy  be 
computed  by  taking  the  first  term  in  a  faylor  expansion 

5  sHipi)  *  ^  x  pa  •  ~  -p((  .. .  (3) 

a  a 
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This  may  be  rewritten  as 

6  <K(p})  "  ±<5<p  <  lp>  |  n  •  L| > 


(4) 


,3  *  ir  and  L  *  I  r  *  p  .  This  may  be  integrated  to 

where  we  have  used  ^  a  tx  a 

give 


4>({p}R(n,4>) )  -  <  (pH  e1<i>n  ‘  Ll*  ” 


(5) 


which  describes  the  effect  of  a  finite  rotation,  R(n,$),  through  an  angle 
4>  about  the  axis  n.  We  have  given  this  derivation  to  emphasize  the  fact 
that  Eq.  (5)  is  a  symbolic  operator  characterization  of  a  generalized 
Taylor  expansion  and  does  not  depend  on  the  nature  of  the  function 
(except  that  we  assume  the  expansion  converges) . 


If  the  rotation  R  corresponds  to  the  rotation  R^  (■  R(n^,82)  followed 
by  the  rotation  R2  =  (R(n2,02))>  then 


<  (p  }R |  tp 


tp)  I 


in-*  L0-  in~*  L0? 

*  1  e  1  I*  >  • 


(6) 


Since  we  wish  to  use  the  symbols  R  both  in  abstract  arguments  and  as 
symbols  for  the  operators  e^n  ^  we  shall  denote  the  product  of  two  rota¬ 
tions  in  such  a  way  that  the  symbols  are  ordered  from  left  to  right, 

R  -  RxR2  . 


We  shall  make  frequent  use  of  the  identity 

Rj-1  R(n,8)  Rx  -  R(nRlt8)  (8) 

which  may  be  proven  by  simple  geometric  arguments  or  by  manipulation  of 
the  operator  expressions  using  the  commutation  relations 

L  x  L  -  il  .  (9) 
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We  shall  make  frequent  use  of  the  irreducible  representations  of  the 
rotation  group  defined  by 

D\r,M  (R("»6))  "  <  L.M'!ei9n'L|L.M  >  .  (10) 

Since  any  rotation  can  be  expressed  as  a  rotation  about  the  z  axis  followed 
by  a  rotation  about  the  y  axis  followed  finally  by  another  rotation  about 
the  z  axis  (Euler  angle  representation)  it  is  only  necessary  to  give  a 
formula  for 

dLM.>M(6)  *  <  L.M' |eieLy|L,M  >  .  (11) 

If  the  usual  phase  conventions  are  adopted  for  the  matrix  elements  of  L, 
then  Eq.  (11)  leads  to  the  formula  given  by  'tgner. 

DERIVATION  OF  CENTRAL  RESULT 

For  a  three-body  sy3tem  interacting  through  two  body  central  forces 
the  time  independent  Schrodinger  equation  is 


P1,P2*P3^  > 


2 

/  dg  V12(q  )  <  P1  +  q,  P2  ~  q,  P3I<|>  >  +  ... 


(12) 


where  the  omitted  terms  involve  and  V^.  The  total  linear  momentum 
is  a  constant  of  the  motion  and  we  shall  work  in  the  frame  where  this 
vanishes.  The  three  momenta  then  form  a  triangle  and  the  wave  function 
depends  on  the  magnitude  of  the  momenta  and  on  the  orientation  of  the 
triangle. 


We  specify  the  orientation  of  the  triangle  by  indicating  the  rotation 

R  which  carries  the  momentum  triangle  from  some  specified  reference  orienta- 

P 

tion  p  to  its  actual  position.  We  shall  introduce  three  standard  reference 
orientations  labeled  by  p  *  1,2,3.  These  are  defined  by  requiring  that  the 
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three  momenta  lie  in  the  X-Y  plane  ano  that  the  directed  normal  is  in  the 
z  direction.  For  P  ■  j  ■  1*2,3  the  momentum  Pj  if  in  the  y  direction. 

These  standard  reference  orientations  are  illustrated  in  Fig.  1.  The  wave 
function  at  the  point  (p)  is  related  to  its  value  at  the  reference  position 
W  by  the  operator  equation 

<  {p}|i|<  >  ■  <  (p}p|Rp|iJ'  >  .  (13) 

The  potential  energy  terms  change  both  the  magnitude  of  the  particle 
momenta  and  the  orientation  of  the  momentum  triangle.  These  effects  e- 
Illustrated  ii.  Fig.  2  for  the  term  The  momentum  triangle  rotate* 

about  the  p^  axis  by  the  angle  $  and  the  magnitude  of  the  momentum  p^  and 
P2  is  changed  to  p^ '  and  P2'  . 

Let  Rp'  be  the  rotation  which  takes  the  momentum  triangle  p^'  P2* 
p^  from  its  reference  position  P  to  its  final  position.  It  is  necessary 
to  find  a  relationship  between  Rp'  and  Rp.  The  rotation  Rp'  can  be  per¬ 
formed  in  three  steps.  First  rotate  the  triangle  p^p^p  *rom  its  reference 
position  P  to  the  position  given  in  Fig. 1(c).  We  denote  this  rotation  by 
R  ’  .  This  is  then  followed  by  the  rotation  R_  which  takes  p,  into  its 
final  position  and  the  plane  of  the  momentum  triangle  into  the  P^P2P3  plane. 
This  is  then  followed  by  a  rotation  through  an  angle  $  about  the  p^  axis. 

Thus 


v  -  rp;3  r3r<p3*^  • 


By  using  Eq.  (8)  this  can  be  rewritten 


V  '  Rp',3  r<*-«  R3  ' 


Similar  reasoning  leads  to  an  expression  for  Rp, 


R  »  R  _  R 
P  P ,  3  3 


(14) 


(15) 


(16) 
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By  eliminating  we  obtain 


R 


P 


t 


Rp;3  R(y,*> 


R 

P 


where  R. 

3»P 


(17) 


It  is  also  necessary  to  change  variable  from  q  to  p1  ' ,  p2',  and  4>. 
The  Jacobian  of  this  transformation  is  most  easily  found  in  two  steps.  We 
introduce  the  symbol  A  for  the  area  of  the  triangle,  so  that 

16AZ  *  W+P2+P3^~Pl+P2+P3^PrP2+P3)(pl+p2~P3)  *  (18) 

Then  2A'/p3,  P3‘P1'/p3,  and  $  form  a  set  of  cylindrical  coordinates.  Thus 

dq  =  (4/p33)  A'dA'd(p1' *P3)  d$  .  (19) 


It  is  then  relatively  easy  to  make  the  additional  transformation  to  obtain 
dq  -  (l/|Pl+p2|)  Pl'  p2*  dpL'  dp2*  d  $  .  (20) 

2 

An  expression  for  q  Is  also  required.  After  a  little  algebra  one  finds 
q2  -  (l/p32)  [(qxp3)2  +  (q*p3)2] 

“  j(4/Pl+P2)2 j  [(A2+A'2  -  2AA’  cos  *)  +  (l/16)(Pl2+p2'2-Pl'2-p22)2]  . 

(21) 


Substituting  this  into  Eq.  (12)  we  obtain 


2  2  2\ 

P1  p2  p3  \ 

2MX  "  2M2  '  2M 3  i  <  {p}p|Rp*  ♦  >  " 


CiPl+PlT)  /pl'dpl'P2'dp2'd*  V12(q2)<fp')oil<D,3  + 


(22) 
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We  have  not  yet  used  the  tact  that  is  an  eigenstate  or  the  angular 

2  2 
momentum  operators  L  and  L  .  If  is  an  eigenfunction  of  L  and  L  with 

z 

eigenvalues  L(L+1)  and  M  it  is  possible  *-o  express  the  matrix  elements  of 
the  rotation  operators  in  terms  of  the  functions  D  . 

The  functions  f  defined  by 

S«P>  (t,M  '  {p,pUL.M  ”  (23) 

depend  only  on  the  invariants  jpjJ»|p2l  an^  I P3  j  •  They  then  satisfy  the 
equation 

(’  2  2  2 

E  1, 

2Mi  2M2  ?M3 

Tk.  (1/ |p1+P2  I )  x  J Pj_'  P2‘  dp^‘  dp2'  <4  V12^2> 

“k'.kMj  *<*•♦>  R3,p>  “  £L,K'(Pl'>P2,-l>3)  +  "•  (24> 

The  relationship  Detween  the  wave  function  at  an  arbitrary  point  {p} 
and  the  f'e  is  given  by 

"  {p}K,M  >  "  <  {p}plRp^L,M  >  "  6(Ep)  h  fL,K(pl»P2»P3)I4,M(Rp)* 

(25) 

It  is  a  simple  exercise,  which  we  leave  to  the  reader,  to  demonstrate  that 
Eqs.  (24)  and  (25)  imply  that  the  wave  function  <  {p)|ip  >  is  independent 
of  what  choice  is  made  for  the  reference  orientation  p  although  the  functions 
fp  do,  of  course,  depend  on  this  choice. 

If  the  reference  orientation  p  is  such  that  the  momentum  triangle  is 

in  the  x-y  plane  then  R  .  is  a  rotation  about  the  z  axis  so  that  D  occur- 

P  *J 

ring  in  Eq.  (24)  is  a  simple  phase  factor  times  d  defined  in  £q.  (11) — that 
is,  we  have  the  usual  Euler  angle  representation  of  the  rotation  group.  In 
particular  2  “  R(z,tr+y),  R2  3  “  R(z,7Tta),  and  R^  R(z,-:f+8). 
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It  should  be  noted  that  Eqs.  (22),  (23)  ari  (24)  are  true  for  any 
number  of  particles  if  we  add  the  needed  additional  kinetic  energy  terms 
on  the  left-hand  sides  and  provide  the  additional  scalar  arguments  for  the 
functions  f.  There  are  however  two  fairly  serious  geometrical  compile 
tions.  The  variables  p^.Pj,  •••  PN  are  no  longer  a  complete  set  of  scalars 
and  an  additional  3N-9  scalars  must  be  introduced.  Some  of  these  additional 
scalars  will  change  due  to  the  action  of  V^2  and  will  be  complicated  func¬ 
tions  of  p^'  ,  p2 '  ,<(>  and  the  initial  set  of  scalars.  In  addition  the  de¬ 
termination  of  the  rotations  R  .  will  involve  rather  difficult  geometrical 

P » 3 

considerations. 


EFFECTS  OF  PARITY,  EXCHANGE,  AND  TIME  REVERSAL  INVARIANCE 

For  a  three-body  system  it  is  fairl,  easy  to  determine  the  effects 
of  a  parity  or  an  exchange  transformation  of  the  functions  f. 

For  a  three-body  sy3tem  the  parity  transformation  is  equivalent  to 
a  rotation  through  180°  about  the  normal,  n,  to  the  momentum  triangle. 

This  is  equivalent  to  first  rotating  the  reference  state  by  180°  about  the 
z  axis  provided  the  reference  state  lies  in  the  x-y  plane.  For  even  (odd) 
parity  states  we  obtain 

1  <  { P }  [  Ip  >  =  i  <  {p}p  jRp  |t|>  > 

■  <  {-p}|4>  >  ■  <  {pip  |  R  ( z ,  7T )  Rp|i|)  >  .  (26) 


This  immediately  leads  to  the  relation 


(PARITY) 
-  (~1)K 


(27) 


If  particles  1  and  2  are  indistinguishable  it  is  desirable  to  know 
the  effect  on  the  f's  of  their  exchange.  For  this  purpose  it  is  convenient 
to  use  the  reference  orientation  illustrated  in  Fig.  1(c)  where  p^  is  along 
the  y  axis.  The  effect  of  an  exchange  is  illustrated  in  Fig.  3  and  it  is 
seen  that  the  directed  normal  to  the  plane  points  in  *-he  -z  direction 
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instead  of  the  ?.  direction.  This  can  be  modified  by  performing  a  rotation 
about  the  y  axis  of  180°.  For  states  even  (!)  or  odd  (-)  under  exchange 
we  find 

’  fL,K^Pl,P2’P3)  "  EK'  f*3  k'  (p2,Pl’P3^  LK *  | R ( y » >  1 LK  >  (28) 


Using 


°K\K  *  V,-K  ("1)L_K 


(29) 


we  find 


*  fL,K(pl’P2’P3)  "  (“1)L"k  fL,-K  (P2»P1»P3)  ‘  (30) 

if  all  three  particles  are  equivalent  it  is  possible  to  use  equations 
analogous  to  Eq.  (30)  for  the  exchange  of  any  pair.  These  results  com¬ 
bined  with 

fL,K  =  EK’  fL,K'  DK'K 

provides  all  the  needed  information  —  unfortunately  in  a  rather  cumbersome 
form.  Attempts  by  the  author  to  define  a  reference  orientation  which 
treats  all  three  particles  symmetrically,  and  which  leads  to  simpler  ex¬ 
pressions  for  the  effects  of  exchange,  have  been  unsuccessful. 


The  final  symmetry  transformation  we  consider  is  the  behavior  under 
time  reversal.  Following  Wigner  we  note  that  solutions  can  be  chosen  such 
that 

<  {-p}i*L,M  =  <_1)M  ''  {P'K,-M  *  *  (31) 


If  we  apply  E;.  (26)  and  observe  that 


(-1) 


M-K 


(32) 
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we  Immediately  obtain 

TIME  REVERSAL 

(f£)*  ■  f"K  03) 

Since  we  have  used  Eq.  (26),  the  reference  orientation  P  must  be  such  that 
the  momentum  triangle  lies  in  the  x-y  plane.  For  bound  state  problems  the 
overall  phase  of  the  wave  function  can  always  be  chosen  so  that  Eq.  (33) 
is  satisfied.  For  scattering  problems  it  is  also  necessary  that  the  inhr no 
geneous  terms  be  chosen  so  as  to  satisfy  Eq.  (33).  (This  may  always  be 
done. ) 
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ON  THE  SLOPE  OF  ENERGY  CURVES 
AS  THEY  ENTER  THE  CONTINUUM 

by 

Peter  J.  Redmond 

ABSTPACT 

Let  the  Hamiltonian  H(X)  have  a  bound  state  with 
binding  energy  E(X)  such  that  the  E(X)  ■+•  0  as  X  -*■  X*. 

It  is  shown  that  dE(X)/dX  ■  0  at  X  ■  X*  for  s-state-. 
and  that  dE(X)/dX  ^  0  for  higher  angular  momentum. 

As  an  application  of  this  theorem  it  is  shown  that 
there  is  a  low  lying  p-wave  resonance  in  the  scatter¬ 
ing  of  an  electron  by  molecular  nydrogcn. 

According  to  the  best  available  calculation1  the  lowest  lying  adia¬ 
batic  potencial  curve  for  the  system  H  +  H  intersects  the  ground  state 
potential  at  a  point  to  the  right  of  the  H^  minimum.  It  has  been  sug¬ 
gested^  that  the  two  curves  should  be  tangential  at  the  point  of  intersec¬ 
tion  whereas  the  numerical  calculations  indicate  that  this  is  not  so. 

Which  of  these  two  possibilities  is  realized  is  important  in  deciding 
whether  it  is  reasonable  to  consider  extending  the  H  +  H  curve  beyond  the 
intersection  point  and  interpreting  it  as  a  potential  resonance.  The  same 
question  arises  for  any  A  +  B  system  and  in  many  other  molecular  problems. 

We  shall  show  that  the  essential  features  of  the  problem  are  present 
in  a  very  much  simpler  situation.  We  consider  the  motion  of  a  particle 
governed  by  a  Hamiltonian  containing  a  potential  which  depends  on  a  pa¬ 
rameter  x»  It  is  assumed  that  the  Hamiltonian  Las  a  bound  state  with  energy 
E(X)  over  some  range  of  values  for  X  and  that  E(X)  +  0  as  X  +  X*.  We 

*  Submitted  to  Phys.  Rev. 

1  H.  S.  Taylor  and  F.  E.  Harris,  J.  Chem.  Phys.  39.,  1012,  1963. 

Private  communication,  B.  F.  Gray. 
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shall  show  that  dE(A)/dA  “  0  at  A  *  A*  for  s  states  and  that,  in  general, 
d(E(A)/dA  i  0  for  higher  angular  momentum. 


Consider  the  Schrodinger  equation  (If  =  1) 


1 

2m  r 


r  + 


m±n  + 

2mr2 


V(r,A)] 


♦(r,A)  =  E(A)  $(r,A)  . 

(1) 


For  simplicity  we  assume  that  V(r,A)  =  0  for  r  >  rQ  although  this  is  not 
essential  to  the  argument.  The  slope  of  the  energy  curve  is  given  by 

r 

■— -  J  $2(r,A)  |y  V(r,A)  r2  dr  .  (2 

o 


We  now  consider  the  behavior  of  the  s-wave  solutions  as  E(A)  •+  0. 

At  distances  greater  than  tq  the  wave  function  is  proportional  to 
e  <r/r  where  *  is  defined  by  E(A)  *  -  <2/2m.  Let  the  solution  to  Eq.  (1) 
for  r  *  rQ  be  denoted  by  x(r,A)  where  we  normalize  x(r»A)  to  unity. 


2,  2  . 

X  (r, A)  r  or 


o 


1  . 


(3) 


Then  the  wave  function  x(r,A)  will  approach  a  definite  well  defined  limit 
as  E(A)  0.  (As  a  simple  example:  for  a  square  well  potential  of  range 
rQ  and  variable  depth  x(r,A)  ->  (2/rQ)1/2  sin  (TTr/2rQ) /r . )  The  normalized 
wave  function  then  is  given  by 


4> (r , A)  =  n  x(r,A)  ,  r  <  r 

o 

-<(r-r  ) 

=  N(ro/r)  x(rQ,A)  e  ,  r  >  rQ  .  (4) 


The  normalization  factor  is  readily  obtained 
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N2  -  (1  +  rQ2  x2  (ro,X)/2<) 


(5) 


and 


N2  ■+  2k/x2  (rQ,X)  rQ2  -*■  0  as  E(X)  -*■  0  . 


2 

The  possib'j-ity  that  the  conclusion  N  *  0  could  be  incorrect  because 

X(r  , X)  ■*  0  can  be  ruled  out  by  noting  that  the  logarithmic  derivative  of 
o 

the  r  times  the  wave  function  nust  be  continuou  •  at  r  ■  rQ.  This  gives 


*<vx) 


*  r  *(r,X) 


r=r 


o 


2k  -*•  0 


(6) 


so  that  drx/dr  ■+  0  at  r  ■  r  .  If  x(r  ,X)  also  approached  zero  then  $(r,X)/N 

o  o 

would  approach  a  function  which  was  identically  zero  everywhere.  This 
contradicts  the  assumption  that  a  bound  state  exists  in  the  vicinity  of 
E  =  0. 


In  order  to  illustrate  the  argument  for  higher  angular  momentum  we 
consider  a  p-state.  The  wave  function  is  then  given  by 

<Kr,X)  **  N  x(r,X)  ,  r  <  rQ 

-< (r-r  ) 

=  N(ry/r)  x(rQ»X)  e  (1  +  1/icr) /(I  +  1/*^)  »  r  >  rQ  • 

(7) 


An  expression  for  the  normalization  factor  is  easily  written  and  one 

2  3  2  -1 

readily  finds  that  as  E(X)  ->-0  N  ■+  (1  +  rQ  x  (rQ,X*))  .  The  argument 

is  readily  extended  to  higher  angular  momentum  and 

Nt2-  (l+r 

Since  the  normalization  factor  approaches  a  finite  non-zero  value  the  slope 
dE/dX  given  by  Eq.  (2)  will  generally  be  non-zero. 


3  2 


(r  ,X*)/(2£ 
o 


-1 


l  >  1 


»  ^  .1 


(8) 
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The  above  simple  treatment  of  a  simple  model  can  be  described  in 
more  general  terms  which  indicate  the  wide  applicability  of  the  results. 
Consider  a  many  body  Hamiltonian  which  depends  on  a  varying  parameter  such 
that  the  binding  energy  of  the  system  approaches  zero.  As  this  limit  is 
approached  one  of  the  particles  will  tend  to  separate  from  the  remaining 
particles.  At  points  where  the  separation  distance,  r,  between  these  two 
components  is  very  large  the  wave  function  is  a  product  of  a  bound  state 
function  of  the  particles  left  behind  multiplied  by  ^“(0,$)/ r  ,  where 

l  is  the  relative  angular  momentum  of  the  two  components.  For  l  =  0  the 
asymptotic  wave  function  is  not  normalizable  at  zero  binding  energy  so  that 
the  probability  that  the  separated  particle  is  in  the  range  of  the  forces 
of  the  bound  state  left  behind  goes  to  zero.  For  higher  angular  moment  w 
the  asymptotic  wave  function  is  normalized  so  that  the  separated  particle 
continues  to  interact  with  the  remaining  particles. 

The  result  is  a  manifestation  of  the  angular  momentum  potential  bar¬ 
rier.  The  non-vanishing  of  dE/dX  for  l  +  0  suggests  that  the  energy  curve 
may  be  extrapolated  into  the  continuum  where  it  would  be  identified  as  a 
resonance. 

As  an  illustration  of  these  ideas  we  consider  the  lowest  potential 
curves  of  the  H2  molecule  and  HH~  system.4  The  general  shape  of  these 
curves  is  illustrated  in  Fig.  1  and  the  question  to  be  decided  is  whether 
or  not  the  curves  are  tangential  at  the  point  of  intersection.  It  is 
easily  seen  that  the  arguments  presented  above  applv  to  the  difference  in 
energy  between  the  two  curves.  The  curve  can  be  loosely  described  as 
being  made  up  of  the  molecular  orbitals  (og)  ou»  At  the  point  of  inter¬ 
section  the  extra  electron  (<0  in  the  HH"  system  is  in  a  p-state  relative 


3  There  are  also  terms  with  lower  inverse  powers  of  r  but  these  are 
multiplied  by  powers  of  K  such  that  the  products  do  not  contribute 
to  the  normalization  Integral. 

4  In  application  it  is  important  that  the  long  range  force  experienced 
by  the  odd  particle  is  dominated  by  the  centrifugal  barrier.  This 
is  satisfied  for  the  e"  -  H2  system. 
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to  the  H2  molecule  so 
there  should  be  a  low 
by  Hr5 


that  the  curves  cross  obliquely.  This  suggests  that 
?ving  p-state  resonance  in  the  scattering  of  electrons 


Figure  1.  Lowest  Lying  Potential  Curves  of  H2  and  HH 


5  Experiment?  1  evidence  for  such  a  resonance,  as  well  as  referenr.es 
to  previous  theoretical  and  experimental  studies,  is  contained  in 
G.  J.  Schulz  and  R.  K.  Asundi,  Phys.  Rev.  Letters  15.  946,  (1965). 
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